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Abstract 

In chapter 1 I begin by discussing the basic ideas of quantum field theory 
(QFT). I provide a review of symmetries in physics and then move on to 
discuss the quark model Chapter 2 is a review of lattice gauge theory with 
particular attention paid to lattice QCD. I begin by introducing lattice QCD, 
I then discuss some of the associated problems. I move on to discuss gauge 
fields on the lattice along with free lattice fermions. I then use this to define 
the lattice QCD action. I conclude this chapter by discussing how to repro- 
duce the correct continuum physics. Chapter 3 discusses the basic numerical 
techniques employed in lattice simulations. I review methods for putting 
particles onto the lattice and conclude with a discussion of how to fit the re- 
sulting data. Chapter 4 reviews symmetries of the QCD Lagrangian, various 
forms of symmetry breaking in physics, the PCAC relation, the Goldberger- 
Treiman relation and the spontaneous breakdown of the axial symmetry. I 
move on to discuss sigma models and finally arrive at a basic chiral pertur- 
bation theory. I present research completed with my supervisor C. Allton 
and collaborators A.W. Thomas, D.B. Leinweber and R. Young in chapters 
5 & 6. This work involves making lattice predictions for the hadronic mass 
spectrum using extrapolation techniques based on the predictions of chiral 
perturbation theory which have been developed by the Adelaide group. 
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Chapter 1 



Introduction 



1.1 Why Quantum Field Theory? 



Although quantum mechanics was a pioneering theory, it was apparent to 
all that it failed on many different levels. The most basic failure of quantum 
mechanics is its inability to account for a relativistic system of particles. In 
such a system the number of particles is not conserved. Dirac knew that 
this inconsistency had to be resolved in order to correctly account for a 
real particle process. In 1927 he published a paper The quantum theory 
of the emission and absorption of radiation which was a first attempt at 
unifying the theory of special relativity with quantum mechanics. It was this 
paper that laid the foundations for a quantum theory of fields, all modern 
theories have their roots based in this. Quantum field theory has proved to be 
an amazingly successful framework for building theories of the fundamental 
forces of nature. Its predictions for the interactions between electrons and 
photons have proved to be correct to one part in 10^. Moreover in the 
form of the standard model, it explains three of the four fundamental forces 
of nature, electromagnetism and the strong and weak nuclear forces. The 
standard model only fails to explain the fourth fundamental force, gravity. 
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1.2 The path integral 

The path integral is a very powerful method of quantisation and is of great use 
in QFTs. Here we review a simple example by considering the Hamiltonian 
for a quantum mechanical particle in one space dimension 

2 

H=^ + V{x) = Ho + V (1.1) 

In the Heisenberg representation we may write the transition amplitude as 

{x\t'\x,t) = (x'|e-'^(*'-*)|a;) (1.2) 

If we use the fact that e""*"* = e°e* and insert a complete set of co-ordinate 
eigenstates, 

dxi\xi){xi\ = 1 (1.3) 



between the exponentials, let T = {t' — t) and At = {ti —t), we then have 
{x',t'\x,t)= /rfxi(x'|e-'^(^-^*)|xi)(a;i|e~^''^*|x) (1.4) 



Dividing T into n equal parts (T = nAt) and inserting (n — 1) states in this 
way we have 

{Xn-i\e-"""K-2} ■ ■ ■ {xi\e-""''\x) (1.5) 

For small At the exponentials can be well approximated using only the first 
term of the Baker-Campbell-Hausdorff formula (eq. 2.19) allowing us to 
rewrite the matrix elements as 

(1.6) 

where we have used the fact that V only depends on space co-ordinates. We 
can calculate the remaining matrix element by introducing a complete set of 
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momentum eigenstates, 

^dp\p){p\ = l (1.7) 
and making use of the fact that 

(.|p> = —e^ (1.8) 

By combining the remaining exponentials and completing the square we are 
left with a simple Gaussian integral, performing this gives 



(1,9) 

Hence our amplitude takes the form 



?t-i 

/ W nn W nn 

J\„-iHT 



{x'\e-''''\x) 



^ \ m ) k=0 



(1.10) 

If we now consider the limit of n ^ cxd we see that the exponent becomes 
the classical action for the path x{t) from x to x'. 



n-l 2 



EMi(^)-''w} 

k=0 

asn^oo £*{|(|)^V(.)}=5 (1,111 



Finally we note that the integrations over the Xk can be interpreted as an 
integration over all possible paths x{t). To describe this we introduce the 
notation 

( ) dxi ... dxn-i ^ const. J~\ dx(t) = Vx (1-12) 

\27rih''At/ 
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Hence we may now write our quantum mechanical amplitude in the path 
integral representation as 



iS 



(x'|e-'^^|x) = y Dxe^ (1.13) 

To make the transition to classical mechanics we simply take the limit h 0. 
To make the transition to a three dimensional theory we simply generalise 
to paths Xi(t) 

Vx-^Y[Y[dxi{t) (1.14) 

t i 

1.3 Quantum Field Theory 

As discussed at the beginning of this chapter quantum field theory is the 
most successful frame work for describing the sub-atomic world. In section 
1.2 we derived the path integral for a simple one dimensional quantum me- 
chanical system. To move to a quantum field theory we must introduce the 
functional integral representation of quantum field theory. Although this can 
be derived rigorously, here I will motivate it by analogy. The key concept is 
to promote the basic variables, Xi{t), of our quantum mechanical example to 
fields, ■ip{x,t). Our rules for the transition are then: 

Xi ip{x,t) 
i —>■ X 
Yl dxi (t) JJ rf^/; (f , t ) = 

t,i t,x 

S = J dtL S = j dSC (1.15) 

Here L is the Lagrange function and C is the Lagrangian density, which 
from here on will be referred to the Lagrangian. The objects of interest in 
quantum field theory are the vacuum expectation values of field operators, 
also known as correlation functions or Green's functions. These Green's 
functions contain all physical information about the system. In analogy with 
our quantum mechanical path integral we can write a representation of the 
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Green's functions in terms of functional integrals: 

(0|O|0) = ^1 V^jOe 



iS 



with Z = J V^pe- (1.16) 
We interpret this as an integration over all classical field configurations. 

1.4 Symmetries 

One of the major advantages of the Lagrangian formalism of QFT is that sym- 
metries of the Lagrangian lead to conserved currents, also known as Noether 
currents. To exemplify this we consider a Lagrangian that is symmetric under 
some given transformation of the fields: 

i/j — > ip + 5ip 
d^^ ^ d,^ + Sid,^) (1.17) 

For a symmetric Lagrangian we have: 

Hence SC = C{il) + d^^j + 5(9^^)) - d^^) = 

where we have Taylor expanded the first term to leading order in Sip. Using 
6{dfj,ip) = d^{ip + 5ip) — d^ip = d^{5ip), the equations of motion for a field"^ 
and the rule for differentiating a product we have: 

= d^J" (1.19) 



^For example see chapter 1 of Quantum Field Theory by Michio Kaku. 
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Here J'^ is the conserved current. With this their is an associated conserved 
charge. To calculate this we integrate the conservation equation over all 
space: 

= / d^xd^J" 

d^xdoJ^ + J d^xdiX 

= d^xJ° + J dSif 

d 

= -^Q + surface term (1.20) 

Assuming that our field [ip) vanishes at infinity, the surface term can be 
neglected. Hence a conserved current leads to a conserved charge. We will 
see the importance of these ideas in chapter 4. 



1.5 The Quark model 
1.5.1 The Eightfold Way 

Oppenheimer once quipped "The Nobel Prize should be given to the physi- 
cist who does not discover a new particle" . He was referring to the seemingly 
endless discovery of new particles that was taking place during the 1960 's. 
Theoretical understanding of elementary particles during this period was a 
mess. Although Yukawa proposed a theory describing the strong interaction, 
it had a coupling constant that was very large and hence perturbation theory 
was unreliable. One important observation was that the existence of reso- 
nances usually indicated the presence of bound states. This lead Sakata [1] to 
postulate that the hadrons^ were composed of states built out of proton (p), 
neutron {n) and lambda (A) particles. Ikeda, Ogawa and Ohnuki took this 
idea further by proposing that these particles transformed in the fundamen- 
tal representation (3) of SU{3). They also stated that mesons could be built 
out of bound states of 3 and 3 [2] . Unfortunately some of their assignments 

^The name hadrons comes from the Greek word hadros meaning strong. 
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were incorrect though. The correct SU (3) assignments were discovered by 
Gell-Mann and Ne'eman. They postulated that baryons and mesons could 
be arranged in what they called the "Eightfold way" [3] . Gell-Mann went on 
to propose (with Zweig) that the SU (3) assignments could be generated by 
introducing new constituent particles called "quarks" which transformed as 
a triplet 3. 



1.5.2 Strangeness 

It had been observed that a new quantum number, in addition to the isospin 
quantum number, was also conserved by strong interactions. This was called 
strangeness, and could be explained in terms of the SU{3) flavour group. This 
group has representations labelled by two numbers, the third component 
of isospin (Is) and a new quantum number called hypercharge (Y). The 
strangeness quantum number and hypercharge can be related to each other 
via the Gell-Mann-Nishijima formula [4], [5]: 



Q 



/ I \ 

-| 

Vo -1/ 



'1.211 



with Y = B + S , where B is the baryon number, S is the strangeness, and 
Q is the charge. 



1.5.3 A global SU{?>) symmetry 

To fit the known hadronic spectrum of particles, it was proposed that mesons 
were formed from a quark and anti-quark, while baryons were formed from 
three quarks. Hence it was expected that mesons and baryons would be 
arranged according to the following tensor decompositions: 

Meson 3 ® 3 = 8 © 1 

Baryon 3®3®3 = 10©8©8©1 
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To see how the bound states are constructed for the mesons we arrange the 
meson matrix according to their quark wave functions: 

/ uu ud us ^ 
M = 3 ® 3 = du dd ds 

\ su sd ss J 

( {2uu — dd — ss)/3 ud us \ 

du {2dd — uu — ss) / 2) ds 

\ su sd {2ss — uu — dd)/3 / 

+ {l/3)l{uu + dd + ss) (1.23) 



Using this we can write the meson matrix for the pseudoscalar mesons as: 



M 



v 



1 



K+ \ 



TT 



;i.24) 



This octet may be represented graphically by plotting isospin against hyper- 
charge. Figure 1.1 depicts this. 



1.5.4 QCD 

After many decades of confusion QCD emerged as the best candidate to 
describe the strong interaction. It has six flavours of quark in the funda- 
mental representation, these can be arranged into three families {u, d), (c, s) 
and (t,b). Leptons can similarly be grouped into three SU(2) doublets in 
electro-weak theory. It is unclear why there should only by three families in 
the standard model. QCD is based on the SU{3) colour symmetry group. 
The eight generators of the group are represented by the Gell-Mann matri- 
ces^ \a, (a = 1, ... ,8). The gauge fields (gluon fields) are denoted by A^. 
We express the gluon field strength tensor as: 

G% = d,Al - d^Al + gfabcAlAl] (1.25) 

■^For example see Quarks Leptons & Gauge Fields by Kerson Huang 
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Figure 1.1: A graphical representation of the pseudoscalar meson octet. J3 
represents the third component of isospin and Y is the hypercharge. 



The quarks are coupled to the gluon fields via the covariant derivative: 

D^ = d^ + igXaAl (1.26) 
Putting this together we have the QCD Lagrangian: 

1 ^ 
Cqcd = —^Gl^G'^^^ + - m;,)z^, (1.27) 

where the Yang-Mills field carries the SU (3) colour force. The gauge group is 
unbroken and hence the force mediators (gluons), are massless. The quarks 
{ipf,h) carry a flavour index (/, h) along with a colour index and Dirac index 
which I have suppressed. 
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1.6 A note on units and notation 

Throughout this document we choose the natural system of units'^ in which 
h = c = 1. We do this to simphfy formulae and calculations. We may move 
back to conventional units via the following: 

h = 6.58 X 10~=^2 [MeV sec] 

he = 1.97x10"^^ [MeVcm] (1.28) 

Another useful conversion factor is: 

he = 197 [MeV fm] (1.29) 

Which we shall employ when setting the scale in our simulations. Throughout 
this document we will frequently employ "slash notation" this is used because 
the product of the Dirac matrices with a four vector occurs so frequently. In 
the Minkowski metric it is defined by: 

7% = 7%o + = ^ (1.30) 



"^This is a system where one unit of velocity is c and one unit of action is h 



Chapter 2 
Lattice QCD 



In this chapter I review some of the fundamentals of lattice QCD. Detailed 
accounts of this subject can be found in [6] & [7]. 

2.1 An introduction to lattice gauge theory 

Quantum Chromodynamics (QCD) is the leading candidate for a theory of 
the strong interaction. Unfortunately perturbation theory fails to reproduce 
nearly all of the low energy features of the hadronic world, an example of 
this would be the spectrum of the low lying hadron states. Perturbation the- 
ory only seems to be effective in the asymptotic region where comparisons 
between theory and experiment can be made. Non perturbative methods 
have proved to be very difficult in Quantum Field Theories (QFTs). One of 
the most powerful and elegant non-perturbative methods is Wilson's Lattice 
Gauge Theory. In principle lattice gauge theory allows us to put QCD on 
a computer and calculate the basic features of the low energy strong inter- 
action spectrum. This approach is only limited by available computational 
power. Monte-Carlo methods have proved very effective in producing pre- 
dictions that roughly match experiment, and with computational power on 
average doubling every eighteen months the discrepancy between theory and 
experiment is ever decreasing. 
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2.2 The price we must pay 



Putting QCD on a lattice comes with a price, 

1. The metric is Euchdean. This means that lattice gauge theory calcu- 
lations are limited to the static properties of QCD. 

2. Lattice gauge theory explicitly breaks continuous and rotational invari- 
ance because space-time is discretised. 

3. Lattice gauge theory is limited by available computational power, so we 
must work with quark masses that are far greater than actual physical 
masses. This also puts constraints on the volume of space-time that 
we can work in. 



Some of these problems can be overcome by taking the continuum limit, 
this is where we let the lattice spacing (a) 0. We must also take an 
infinite volume limit. However the limitation of computational power means 
that currently lattice sizes are modest. 



2.3 The path integral on a lattice 

We begin by making a Wick rotation. Put simply if {x'^ , , , x^) are co- 
ordinates in Minkowski space-time (with x^ being the time coordinate) then 
we set: 



x^ = ix"" (2.1) 

This can be thought of as a rotation in the complex time plane and gives 
us a imaginary value for our time coordinate. The new set of coordinates 
{x^ , x'^ , x^ , x^) now have a Euclidean metric. The main benefit of doing this 
is that the action (S) is now a real positive quantity and our phase factor (see 
eq. 1.16 becomes a real weighting and so can be interpreted as a probabihty. 
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2.4 Space-time discretisation 

Lattice QCD relies on a discrete space-time. Discretising space-time removes 
the infinite number of degrees of freedom available to the fields and replaces 
them with a finite number. This allows the path integral (sec. 1.2) to be given 
an exact definition. We formulate our theory on a hyper-cubic lattice using 
the Euclidean coordinates ( ). Our lattice is typically defined by 

all of the points that obey: 











a 


where 


n* 




G 


Z 


with 





< 


rt < 


L for z = 1,2, 3 


and 





< 


< 


T 



The spacing between lattice sites is known as the lattice spacing, (a), and has 
dimensions of length. L is defined to be the length of the lattice measured 
in lattice units and is a dimensionless number. We apply periodic boundary 
conditions to the spatial dimensions and anti-periodic boundary conditions 
to the time dimension. This ensures Fermi-Dirac statistics. In doing so the 
momentum space is discretised, we have: 

= 

Here the have the same constraints as before. The beauty of discretising 
space-time is that there is now a maximum allowable momentum. This means 
that lattice gauge theory has an ultra-violet cutoff and hence gauge theories 
on the lattice are naturally regularised. 

2.5 Gauge fields on the lattice 

We begin by defining a link between two neighbouring sites on our hyper- 
cubic lattice. We allow each link to have a dynamical degree of freedom which 
we denote by U{n, n + fi) = Uij where fi is the unit vector in the fj, direction. 
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The dynamical degree of freedom Uij belongs to the compact group^ G, for 
example: 



Z2 : 
f/(l): 

SU{N): U,, 



±1 



(2.2) 



: X matrix with det U = 1 k W = U^^ 
We note that the link has an orientation: 

4 = u^' = u,. 



(2.3) 



Hence taking the inverse reverses its direction. 

At lattice point n we define the simplest closed path on the lattice, this 
is the plaquette. It is illustrated in fig 2.1. 



n + u 



n + fi + 




n n + fi 

Figure 2.1: The plaquette is the simplest closed path on the lattice. 



Mathematically we have 

Up = U{n, n + (i)U{n + fi,n + fl + i>)U{n + fl + i>,n + i>)U{n + 0, n) (2.4) 

Using this we may now define the Wilson action for the gauge fields. This is 
just the sum over all distinct plaquettes, P. 



Sg = l3^{l-ReTiUp) 

p 

^The local gauge symmetry group for QCD is 5'[/(3) 



(2.5) 



2.6. FREE LATTICE FERMIONS 
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For lattice gauge theory to predict properties of QFTs we must perform 
functional integrals (eq. 1.16) for example; 

(O) = ^JVUOe-^o 
where Z = J VU e^^« 

Here VU is analogous to VA in eq 1.16. VU is called the Haar measure, it 
is defined as the product over all links: 

Pt/ = JJ dUij (2.6) 

The Haar measure is a way to assign an invariant volume to subsets of locally 
compact topological groups. It has the following properties 

JdUf{U) = JdUf{QU) = JdUfiUQ) withfiGG 

=> left-right invariant. 

J dU = 1 ^ Normalisable. 

J dU U = => Expectation value of 

a gauge non-invariant 
object is zero. 

One important point to note is that on the lattice the volume of the gauge 
group is unity hence no gauge fixing is needed. 



2.6 Free lattice fermions 

I will now discuss the discretisation of the fermion fields. As we shall see 
this must be performed carefully to avoid the dreaded "fermion doubling 
problem" . We begin by considering a naive discretisation of the free fermion 
action: 

Sf=[ d^x ^(x)(7^9^ + m)iP{x) (2.7) 
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We note that the four dimensional integral may be represented as a sum as 
follows, 

jd^x-^a^Y^ (2-8) 
A symmetric difference approximation for dfj,ip{x) is: 

dM^) = "^^^^"/"-^ (2.9) 

Using eq 2.8 and substituting eq 2.9 into the free fermion action (eq 2.7) 
leads us to the lattice action for free fermions. 

3 4 

Sf = J2[Y^ ^nT'^l^n+A - V'n-A) + ma^V'nV'n] (2.10) 
n fj,=l 

To calculate the lattice propagator we use a Fourier transform to move to 

momentum space. We find the lattice propagator in momentum space is 
given by: 

Propagator^^ ~ —sm{ak^) +m (2-11) 

This propagator has bad behaviour as we take the continuum limit (a — > 
0). As we expect, the lattice propagator has a node at A; = 0, but it also 
has a node at the edge of the Brillouin zone {k^ = ^) for each /i. Hence 
the naive discretisation prescription has an unphysical doubling problem for 
each space-time dimension. Wilson proposed a convenient solution to this 
problem. He suggested that the lattice fermion action should be modified by 
hand. We may do this as long as the correct continuum limit is obtained. 
We add the following Wilson term to our previous naive action: 

^Ipni'ipn+il + 'ipn-fL-'^'ipn) (2.12) 

2a 



Calculating the momentum space contribution by Fourier transforming the 
fermion fields for this term and adding it to our previous naive action 
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gives^: 



(W^^^ ^ ^E^'^^+^-E^v — ^^^^ ^^-^^^ 



The new cosine term preserves the minimum at /c = but ehminates the 
unwanted minimum at the edge of the Brillouin zone. This solution to the 
doubhng problem does come with a price, and that is the Wilson term breaks 
chiral symmetry at finite lattice spacing. 



2.7 The lattice QCD action 

The lattice QCD Lagrangian contains the following fields: 



^lj(n) = [^cq]('^) Quark fields having fiavor: / = !,..., Nf 

Colour: c = 1, 2, 3 & Dirac index: a = 1, ... ,4 
U{n, n + ft) G SU (3) link variables: /i = 1, . . . , 4 

It is constructed to be invariant under an SU{3) gauge transformation. The 
gauge transformations for the fields dbTQ clS follows 

U'{n,n + (j) = Vt{n) U{n,n + (1) Vl\n + (x) 

^'(n) = ij{n) n^{n) (2.14) 
where the fl{n) are SU{3) matrices. The Wilson action for QCD is 

SQCD = Sg[U]+Sg[U,4^,ij] (2.15) 

The gauge action {Sg[U]) is given by equation 2.5 and the quark action 
{Sg[U,ip,i>]) is given by 

S,[U,^,^lJ]=^|J^K^y^|Jy (2.16) 



^Where we now use the hash (#) notation to represent lattice quantities 
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where K2;j/[[/] is the quark matrix. The x,y indices represent space-time, 
colour, spin and flavour. The quark matrix is given by equations 2.10 & 
2.12, and introducing a gauge interaction. 



K 



[U] = Sxy -kY, [Sx,y-f.{r - 7M)t^x,M + S^-i^Ar + 7M)f^i J (2-17) 



Wilson's choice for r is one. and in this case there is no species doubling. 
The spinor indicies are carried by the gamma matrices, the colour indices by 
the link variables and there is a Kronecker delta in flavour space, all of which 
are suppressed. The hopping parameter (k) is related to the free quark mass 
by 

K= — ^ (2.18) 

2m + 8 ^ ^ 



2.8 The continuum limit 

For our formalism of lattice gauge theory to be correct it must reproduce 
the correct continuum physics when we take the continuum limit. Here we 
briefly outline this for the gauge part of the action {Sg[U]). The first step in 
doing this for G = SU (N) is to use the fact that a unitary matrix may be 
written as the exponential of an imaginary matrix: 

U^{n) = exp (za^^f v4-(n)) 

Where g = Coupling constant 

A" = generators of the gauge group 
^/^(^) = The gauge fields 

We then use the Baker- Campbell-Hausdoff formula to rewrite our plaquatte. 
Up (eq 2.4), as a single exponential. 

^A^B ^ ^A+BmA,B] + ... ^2.19) 
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We then use the fact that in continuum 

[A^{x),A^{x + i>)] = [A^{x),A^{x)] 

[A,{x),A,{x + jl)] = [A,{x),A,{x)]=0 (2.20) 

to identify the combined exponent with the Yang-Mills field strength tensor: 

F^,{x) = d^A,{x) - d,A^{x) + ig[A^{x), A,{x)] (2.21) 

We then take the trace of the plaquette by expanding the exponential. Next 
we define the Lie algebra for G 

[A'^,A^] = 2i/°^^A^ 
tr(A° A'^) = 2(5°^ 

(2.22) 

and use it to show 

tr[F,^{xr] = If^^F^-^- (2.23) 

We then use tr{Up) and equation 2.8 to rewrite the lattice action. Finally, 
it can be shown that the correct continuum physics expression is reproduced 
as we let a ^ 0. 



2.9 Setting the scale 

Throughout this chapter our discussion of lattice QCD has been in lattice 
units, i.e. we rescale the fields and masses by appropriate powers of the 
lattice spacing to render them dimensionless. This is of no real use if we 
wish to make physical predictions that we may compare with experimental 
values. To be able to do this we must give our lattice predictions their 
correct dimensions. This is called setting the scale. The continuum value of 
an observable (Ocont), is given by 
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where is our lattice observable and is the energy dimension of Ocont- 

Massless lattice QCD contains one free parameter (/?), so we use one ob- 
servable to determine the lattice spacing (a). We may then make physical 
predictions based on our lattice simulations. In massive lattice QCD addi- 
tional observables are needed to set the quark mass. More generally addi- 
tional parameters are needed as more parameters are introduced. 



Chapter 3 



Numerical Methods 



3.1 The effective gauge action 

To compute observables in QCD we define the expectation value of an arbi- 
trary operator O as: 

{O) = — I VUVipV^p O e-^s-^" 



Z . 

Z = I VUVtpVij e"^«~^' (3.1) 



Grassmann variables cannot be modelled on a computer. Hence we cannot 
use any computational method which involves an action that contains Grass- 
mann variables. We can however analytically integrate out the fermion fields 
from the functional integral (Z) . 

{O) = ^ JvU O e^^« det K 
Z = [ VUe-^^deiK (3.2) 



Doing this leaves us with a functional integral over the gauge fields which 
may be expressed as integrals over real numbers. These can be handled on 
a computer with relative ease. 

We are now ready to introduce an effective action, 5'e//, by making use 
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of the following identity: 

detX = ei°^^*^ = e^'-''^^ (3.3) 
So that we now write the integrand of equations 3.2 as: 

S^jf = Sg- ln(det K) = Sg- Tr In K (3.4) 



3.2 The quark propagator 

The quark propagator is a simple example of integrating out the fermion 
fields leaving us with an object that can be calculated on the lattice. The 
quark propagator is defined as the expectation value of the product of a 
and a ip field. 

ii^xi^y) = ^ J V^.VijyVUe^'^^-^'^i^Jy (3.5) 

As before (sec 2.7) the x,y indices represent space-time, colour, spin and 
flavour. The 5*^ & Sg represent the gluonic and quark parts of the lattice 
QCD action respectively. We note that the propagator is not gauge invariant 
because we may apply independent gauge transformations at each x, y point. 
Performing the fermion integration yields: 

(^x^.) = ^J VUe''^ffK-y' (3.6) 

We may now define the quark propagator, G = K"^, for a given gauge 
configuration U. 

Gui^^ f^)^(/3,2,s/), (7,3,2) [U] = ^aT^ia^z (3.7) 

The Greek indices represent spin, the numbers represent colour and the ro- 
man indices represent space-time co-ordinates. The flavour dependence of 
G is a Kronecker delta and so is suppressed. Repeated indices should be 
summed over. This equation (eq. 3.7) can be solved using a matrix inversion 
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algorithm such as the conjugate gradient method. 

3.3 Monte Carlo simulations 

Lattice gauge theory made a great leap forward when the Monte Carlo 
method was introduced. This is because a naive calculation of the path 
integral is prohibitive, since the sum contains a massive number of terms. To 
exemplify this consider the simplest group that we can define on the lattice, 
Z2 with elements ±1. If our lattice had 8^ sites then the path integral sum 
would contain the following number of terms: 

22" ^ 2^6384 ^ ^05460 (3 

Since the number of links is 4 x 8^ = 2^^. 

The Monte Carlo method is an example of importance sampling. It applies 
certain approximations to the path integral which alleviates this problem. 
Normally the path integral sums over an enormous amount of configurations 
that make an insignificant contribution to the integral. If we could ignore 
those configurations and only sum over the ones where the action is near its 
minimum, then our calculation would be much quicker. The Monte Carlo 
method does exactly this. We define a set of initial values for each link on the 
lattice (Si). Then the Monte Carlo method tells us to generate a sequence 
of configurations S2, S3 . . . such that when statistical equilibrium is reached 
the probability of encountering a particular configuration Sj in the sequence 
is proportional to the corresponding Boltzmann weight, Wi = e~'^t^'l The 
smaller set of configurations that we now use {Sj} are those that are near 
minimum action and hence contribute most to the path integral. 

3.4 The Metropolis algorithm 

A common method for generating the sequence of configurations {Sj} is the 
Metropolis algorithm [6]. Consider generating a new configuration a[ from 
the configuration ai by updating a single link using some random process. 



24 



CHAPTER 3. NUMERICAL METHODS 



It is possible to calculate the change in the action via equation 3.9. 



A random number r is now chosen between and 1. If e~ > r then the 
configuration is accepted, if not it is rejected. If AS' is negative then the 
change is always accepted because e"^"^ > 1. If however we only accepted 
negative values of AS* then the action would be constantly decreasing and 
hence would tend towards the classical equations of motion. Of course this 
is to be avoided because it neglects all quantum corrections. 

By choosing a random number (r), we are actually allowing for positive 
AS, hence the action may increase as we change from ai to a[. This al- 
lows for quantum fiuctuations around the classical equations of motion. The 
algorithm progresses by moving to the next lattice site and changing it in 
some random way. Hence another configuration is generated, we test this 
to see if it meets the proper criteria and move on. In this way we sweep 
through the entire lattice successively making small changes. After many 
sweeps through the lattice we begin to reach thermal equilibrium, this yields 
the set of link variables Si. The process is then repeated and the second set 
of link variables S2 is obtained, and so on. Slowly a set of configurations is 
built up {Si}. The effect of the algorithm is that the new configuration a[ 
is accepted with the conditional probability of e~^'^. 

3.5 The quenched approximation 

As seen in section 3.2, S^jf is the correct action for lattice QCD. Unfor- 
tunately the second term in equation 3.4 makes the action non-local. This 
means that generating a sequence of configurations (sec. 3.3 & 3.4) is far 
more computationally demanding than generating a sequence of configura- 
tions for an action that is local. This is because we have to calculate the 
determinant of the quark matrix in the non-local case. By replacing the de- 
terminant in equation 3.4 with a constant that is independent of the gauge 
fields we have a modified the action which is now local. This is called the 




(3.9) 
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quenched approximation. In practice we set det = 1 in equation 3.4 which 
modifies the effective action so that S^jf = Sg. This corresponds to setting 
the hopping parameter (sec. 2.7) to zero. This leaves us with infinitely heavy 
quarks which do not contribute to the effective action. This effect is coun- 
tered by adjusting the remaining parameters of the theory. Although this 
seems like a very crude thing to do, it works surprisingly well with most of 
the essential features of QCD remaining. Quenched results of calculations 
of the light hadron spectrum are within 10% of experimental results. 

3.6 Hadron correlators 

Correlation functions are used to measure many physical observables on the 



Here 0{x) is an interpolating operator. Any gauge invariant combination of 
fermion fields and gauge links can be used as a interpolating operator. As 
above this multi-dimensional integral is well approximated by the Monte- 
Carlo method (sec. 3.3) allowing us to represent the correlation function as 
the average of the operator 0(a;)0^(y) evaluated on each of the independent 
configurations {Si}: 



Where M.(T,i) is the value of the operator 0{x)0^{y) calculated using the 
configuration Sj, Nc is the number of configurations and i represents the i^^ 
configuration in the sequence of configurations S. 



lattice. 




VUO{x)0^{y)e-^^" 



(3.10) 




(3.11) 



3.6.1 Mesons 



To place mesons on the lattice we use interpolation operators, 0{x), which 
are used to create and annihilate mesons. An operator for the meson state 
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\A4) must satisfy: 

{0\O{x)\M) ^ 

{0\O{x)\M') = (3.12) 

Where is any state that is hghter than state This ensures that 

the operator has a non-zero overlap with the state that it is intended to 
represent and has no overlap with any lighter states. For the operator to be 
gauge invariant we require: 

= 5>r^^f ^ (3.13) 

With a & /3 being the Dirac index, a hh the colour index, f ^ g the flavour 
index and x is the space-time co-ordinates. Since J dQ^^x = it follows that 
the expectation value of a meson propagator is zero unless x = y. To define 
a particular type of meson (pseudoscalar or vector) we pick an interpolating 
operator with the same quantum numbers as that type of meson. It is conve- 
nient to represent mesons in written form using the following combination of 
quantum numbers, J^*^. Table 3.1 gives a description of these numbers and 
their corresponding formulae. Table 3.2 gives examples of the J^*" numbers 



Quantum number Description Formula 



J Total angular momentum J = (l + s,l + s — 1, . . . ,\l — s\) 

P Parity number P = (—1)''^^ 

C Charge conjugation C = (—1)'+'^ 



Table 3.1: This table contains a description of the J^^ quantum numbers. 
The orbital angular momentum eigenvalue is represented by / and the spin 
eigenvalue is represented by s. 

for some well known pseudoscalar and vector mesons. As an example we 
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Meson J 



TT 


0-+ 


P 


1 — 


a 


0++ 


bi 


1+- 


ai 


1++ 


a2 


2++ 



Table 3.2: Some of the well known mesons and their J^'-^ quantum numbers. 

consider the following interpolation operator^: 

Oix,t)=^f{x,th,^g{x,t) (3.14) 

This operator (Eq. 3.14) has the same J^*" quantum numbers as the tt- 
meson (pion). It is a local operator because the quark and anti-quark fields 
are at the same site (x, t). It has been shown that for some mesons such as 
the p-meson an interpolation operator that also includes some form of wave 
function between the quark and anti-quark fields gives a better signal or over- 
lap. This technique is known as smearing. Another technique involves the 
addition of bent paths to the original paths between quark and anti-quark, 
this is known as fuzzing. We can combine these interpolation operators to 
form mesons on the lattice. Substituting these operators into equation 3.10 
represents a physical process in which a pion is created at the space-time 
point (0,0), from there it propagates in space-time to the point (x, t) where 
it is destroyed. This is pictorially illustrated in figure 3.1. 

The correlator can give us lots of information about the particle that we 
wish to study. As an example we will consider the local pion correlator. We 
begin by inserting our interpolation operator (Eq. 3.14) into the correlation 



^The index on the quark fields is a flavor index. 
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Space 




0(a;,t)+ 



Figure 3.1: A pictorial representation of a meson correlator on the lattice. 
The meson is created at position and time (0,0) and is destroyed at {x,t). 



function (Eq. 3.10). We have: 

a(a;,t;0,0) = (0|T^/(x, t)75V^3(x, t)^/(0, 0)75^,(0, 0)|0) (3.15) 

Where we have set {y, t) = (0, 0) and T represents time ordering. Next we 
perform a Wick contraction, 

a(x,0;0,0) = ^ JvUe-'^ffTiiGfiO,xh,Ggix,Oh,) (3.16) 

Here the sea quarks are contained in Seff and the valence quarks are con- 
tained in the quark propagators Gf{0,x) & Gg{x,0) (sec. 3.2). This is an 
example of a point to point correlator. As such its momentum is undefined. 
We may represent this integral using the Monte-Carlo approximation. If we 
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treat quark flavors g and / as degenerate we have: 

1 

C^{x,0) = — ^Tr(G(0,a;)75G(a;,0)75) (3.17) 



We will now use this to study the time slice correlator and show how to 
extract an estimate for the effective mass of the pion. We begin by performing 
a Fourier transform over the three spatial dimensions, x, (i.e. over a time 
slice). This fixes the momentum. 

a(t,p) = ^a(x,t;0,0)e*^-^ (3.18) 

X 

Inserting a complete set of energy eigenstates that have the norm: 

I Qn) {n, qn\ = 1 (3.19) 

gives: 

CAt,p) = I ^$^(0|O(x,t)|n,g„)(n,g„|Ot(0,0)|0)e'^-^ (3.20) 



x,n 



We now make use of the fact that a quantum operator O in the Heisen- 
berg picture is time dependent 0{x) = C)(o)e~*-^"*"*'^"''^ and the fact that 
{A\0\B) = {{B\O^A))* inserting this in Eq. 3.20 gives: 

CAp,t) = J ^5^K0|O(0,0)|n,g„)pe-^-*e-(''"-^)- (3.21) 

x,n 

We may now simplify this by noting g-«{gn-p).a: _ §(^q^_p^ g^j^^ J (i^q^§(^q^_ 
p) implies p = qn- We also rotate to Euclidean space it t. This leads to 



Since we are in Euclidean space the correlator is exponentially damped in 
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time so as t oo the lightest meson state dominates (n = 1). 

To extract a mass prediction we work in the rest frame (i.e. we set p = 0). 
We define the effective mass as: 



Plotting this against time allows us to gain a lattice estimate for the pion 
mass by looking for a plateau in the data. To convert this into a physical 
value we must set the scale (sec. 2.9). In practice the hadron mass is obtained 
by fitting c(0, t) to an exponential over a time range where the ground state 
dominates. 

3.6.2 Baryons 

We place baryonic particles on to the lattice in much the same way as we did 
mesons. To define gauge invariant baryon fields we require: 



We can prove that this is gauge invariant by suppressing all but the colour 
indices, and applying a gauge transformation: 



Gauge invariance follows from eabc^a'^l'^c' ~ det(fi) = 1. 

As for the mesonic case, we construct interpolation operators with the 
quantum numbers of the particle that we wish to study. An example for the 
A"''^ would be: 




(3.23) 




(3.24) 



(3.25) 




(3.26) 



With C being the charge conjugation matrix ip^'^^ = —(C'^^)'^. This is an 
antisymmetric unitary 4x4 matrix that relates to its transpose. Explicitly 
we set = ijj^ = ijj^ = u with u being an up quark. As before (sec. 



3.7. FITTING THE DATA 



31 



3.6.1) we may use these interpolation functions to form baryon correlation 
functions and using similar techniques we can extract physical information 
about baryonic states. 



3.7 Fitting the data 

In any lattice calculation a set of observables Oi are calculated on a number of 
configurations. To extract physical observables we fit the data to a particular 
function or model. This is done by the minimisation of the value. In our 
study we work with uncorrelated data. We define ^-s: 



(3.27) 



Where Xi = independent variable 

Hi = dependent variable 

/(xj, {aj}) = fitting function / model 

N = number of data points 

CTj = uncertainty of the i*'' data point 

The is minimised by the set of independent variables that satisfy equation 
3.28 2 

^ = (3.28) 

If the function /(xj) is linear then equation 3.28 may be solved analytically, 
otherwise we must minimise numerically. In our study we quote values of 
reduced also known as per degree of freedom (x^/rf.o./). This can be 
defined as: 

xVd.o.f = TT^Ef^^^^l' (3.29) 

1=1 

Where m = number of fit parameters 



An acceptable fit should have a /d.o.f of about one. 
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3.8 Statistical errors 

We can associate an error with the fit parameters known as the statistical 
error. In the hmit of an infinite number of gauge configurations we expect 
this to go to zero. Typically a lattice simulation has a small number of config- 
urations (9 (few hundred) this is because configurations are computationally 
expensive to generate, especially dynamical ones (sec. 3.5). Consequently 
it is important to quantify any error due to this. To do this we use boot- 
strap error analysis [8]. We begin by creating a bootstrap sub-ensemble by 
selecting n configurations at random. Next the fitting procedure is applied 
to this sub-ensemble, in the same way it is for the true sample. This gives 
a new estimate for the fit parameters. We repeat this procedure riboot times. 
Doing this generates a bootstrap distribution for each of the fit parameters. 
Representing this as a histogram allows the central 67% region to be defined 
and hence the upper and lower error bars. 



Chapter 4 



An introduction to Chiral 
Perturbation Theory 

In this section we will review the symmetries of the QCD Lagrangian. We 
will then show how these symmetries can be exploited so that we can build 
an effective Lagrangian that can be used to describe the low energy dynamics 
of QCD. 



4.1 Symmetries of the QCD Lagrangian 

Along with the U{1) symmetries that give rise to the conservation of baryon 
number, charge and strangeness (eq. 4.1), there exist two other interesting 
symmetries of the QCD Lagrangian, namely the Vector and Axial- Vector 
symmetries. 

ip — >■ e~'^°'^ip B is the baryon number of ifj 
ip — > e~'^^^ip Q is the charge of ip 

—>■ e~^'^^ip S is the strangeness of ip 

where a, /5 & e are arbitrary real numbers (4-1) 
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4.1.1 The Vector symmetry 



ip ^/;exp( i-r^r ) (4.3) 



We will begin by considering the case of massless QCD. Here we need only 
consider a Lagrangian of the form: 

C-massless = ^JP'ip (4.2) 

I have suppressed all indices and ignored the gauge fields as they will not be 
influenced by the transformations that we will consider. 

The vector transformation (Ay) which belongs to the group SU{2)v is 
defined by, 

^ ^ exp^^-z^r'^r^^/' 
2' 

Where t"- are the Pauli iso-spin matrices. It is immediately obvious that our 
massless Lagrangian (eq. 4.2) will be invariant under this transformation 
since the transformation, Ay, has no space-time dependency. Hence our 
Lagrangian is SU{2) flavour (iso-spin) invariant. 

As was demonstrated in section 1.4, symmetries of the Lagrangian lead 
to conserved currents. They can be calculated using equation 1.19. We find, 

= yV^ (4.4) 

This is the vector current, its associated conserved charge is the isospin 
charge. 

We will now consider adding a mass term. The up and down quarks have 
masses around 10 MeV. The next lightest quark is the strange quark and 
this has a mass of the order of 100 MeV. The lightest hadron is the pion 
with a mass of about 140 MeV. So in the low energy limit of QCD only the 
lightest quarks (up and down) need be considered. To do this we define a 
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mass matrix {mud)-, and write the fermion {ip) fields as an isospin doublet: 





' rriu 










. 






V d ) 







(4.5) 



We use this to write down a Lagrangian involving just the up and down 
quarks: 

Cud = qiJPq - -ji'TTiu + md)qq - -(m„ - md)qT^q (4.6) 

Now because of the Pauli matrix in the last term the Lagrangian is only 
invariant if we assume that the up and down quarks are degenerate in mass 
(m„ = nid). In this case the symmetry leads to three conserved currents 
corresponding to the three generators of SU{2) (the Pauli isospin matrices). 
The corresponding isospin charge operators obey the SU{2) relations: 

= ieijkh (4.7) 

This is exactly the same case as for quantum mechanical spin and so we know 
that the eigenstates and eigenvalues must behave in the same way, hence: 

= /(/+l)|/,/3) 

h\I,h) = h\I,h) (4.8) 



This is our first glimpse at an effective Lagrangian., this is a Lagrangian that 
describes physics in terms of experimental (hadronic) degrees of freedom 
rather than fundamental ones (quarks). However we know that in nature 
isospin invariance is broken due to a finite difference between the up and 
down quark masses. But as we shall see in the following section, if the 
breaking is small compared to the relevant energy scale of the theory, the 
symmetry may be treated as an approximate one. 
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4.1.2 The Axial Symmetry 

As before we begin by considering the massless case. The axial-vector trans- 
formation (A^) is defined by: 

^ ^ exp|^i^75r"r^?/' 

^ V;exp|^z^75r"r^ (4.9) 

This time it is not immediately obvious that the massless Lagrangian (eq. 
4.2) is invariant under this transformation. We must pay special attention 
to the derivative part of the Lagrangian. Since the transformation has no 
spatial dependency we can move the exponential through it. 



^ exp (^TV'^r exp (^7^"^'^) d^i) 



(4.10) 



To deal with the problem of moving the exponential through we Taylor 
expand the exponential and re-express it in terms of trigonometric functions: 

exp(-7Vr)7'' = 7^cos( -r"r ) +i7Vsin( -r"r 



2' /' ' V2 / V2 

i 
2 



7'^exp( --7Vr ) (4.11) 



Where we have made use of the anti-commutation relation {7^, 7^} = 0. 
Hence the massless Lagrangian is invariant under the axial transformation 
(Ayi). Again using equation 1.19 we can calculate the associated axial current, 
we find: 



1 - 



^1 = -^i^l^rr^i' (4.12) 

We now consider the case of adding an arbitrary mass term (m). As in the 
massless case the covariant derivative is symmetric under the axial transfor- 
mation, but the mass term is not. 

bL = -m(^^exp(^^7V'^r) exp(^^7V"r)v^^ 



4.1. SYMMETRIES OF THE QCD LAGRANGIAN 



37 



= -m(^^/;)exp(i7V°r) (4.13) 

As with the vector symmetry if the symmetry breaking term is small com- 
pared to the relevant energy scale of the theory then we may regard the 
symmetry as an approximate one^. 



4.1.3 Chiral symmetry 

Chiral symmetry is often referred to by its group structure SU {Nf)A ® SU {Nf)v 
where the subscript A(V) refers to the axial(vector) symmetry. We have con- 
sidered the case of Nf = 2 but the case with A'^^ = 3 is equally valid. We 
begin our discussion of chiral symmetry by introducing the idea of chirality. 
We define operators to project out the left and right handed components of 
the isospin doublet introduced in section 4.1.1. 

= ^(1 + 7^) 

(4.14) 

Applying these projection operators to the isospin doublet gives: 

qL = ^Lq qR = ^Rq (4.15) 

where q = qL + qR- We can use these chirality states to rewrite the massless 
Lagrangian (eq. 4.2). 

^massless = qL0qL + qR0qR (4.16) 



^Consider a circle in 2D. This is invariant under rotations in the plane. Now imagine 
a flat on the circle. This breaks the rotational invariance, but if the flat is small the circle 
will still look very similar under a rotation. Hence the symmetry is approximate. 
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This Lagrangian is invariant under the independent symmetry transforma- 
tions: 

L" : ql^ exp{iT^a'')qL 

i?'^: g« ^ exp(zr'^/3'^)g^ (4.17) 

Again we apply Noether's theorem (using equation 1.19) to find the associ- 
ated conserved chiral currents: 

R'^" = QRl^lryn (4.18) 

Taking hnear combinations of the chiral currents allows us to express the 
conserved currents in vector and axial-vector form: 

a; = R>^'- - L^'- = q^^^^^q (4.19) 
These are the conserved currents from sections 4.1.1 & 4.1.2 respectively. 



4.2 The chiral transformation of mesons 

Before we construct a chirally invariant model we must first investigate the 
chiral transformation of mesons. We first define quark operators that have 
the same quantum numbers as those mesons that we wish to study. 

Pion: vr" = iipr'^'j^ip 

Sigma: a = %hth 

^ ^7 (4.20) 

Rho: pJJ = ipT^-'yij^i/j 

ai: = ^r'^7/.75^ 



We note that the Sigma particle is not observed in the mesonic spectrum. 
Here we define it to be a particle that carries the quantum numbers of the 
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vacuum. We now individually apply the vector (eq. 4.3) and axial-vector 
(eq. 4.9) transformations, for infinitesimal rotations^ to the pion and rho 
mesons. 

Vector transformation: 



vr" + r X vr" (4.21) 
pa ^ p^ + rxp;: (4.22) 

where, for the pion, we have made use of the commutation relations for the 
Pauli matrices [r°, r^] = 2ieabc'^'^- These transformations are just isospin 
rotations by an amount 6*". 
Axial transformation: 



TT^ + e^a (4.23) 
^ p;: + rx<,^ (4.24) 

this time we have made use of the anti-commutation relations for the Pauli 
matrices {r°,r^} = 25"^ and the fact (7^)^ = 1. These results indicate that 
the axial symmetry is a symmetry of the QCD Lagrangian, and that particles 
rotated into each other should have the same eigenvalues. This is clearly not 
the case since the rho and ai have very different masses. This cannot be 
accounted for by the explicit symmetry breaking due to the mass of the light 
quarks because they have a small mass. This problem is resolved by the 
introduction of the spontaneous breakdown of the axial symmetry. 



4.3 Symmetry breaking in physics 

We will be concerned with two types of symmetry breaking, explicit and 
spontaneous symmetry breaking. 

^If we consider infinitesimal rotations we may Taylor expand the exponential to leading 
order in 9. 
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4.3.1 Explicit symmetry breaking 

Explicit symmetry breaking occurs when a Lagrangian has a symmetry that 
is broken by the addition of some term. Section 4.1.2 exemphfies this, the 
massless Lagrangian (eq. 4.2) is invariant under the axial transforms (eq. 
4.9). But we see that the addition of a mass term (eq. 4.13) breaks this 
symmetry. 



4.3.2 Spontaneous symmetry breaking 

A symmetry is said to be spontaneously broken if the Lagrangian of a system 
possesses a symmetry which its ground state does not. An example this is 
the spontaneous breakdown of a rotational symmetry in a ferro magnet. The 
Hamiltonian for such a system has the form: 



where (ft represents the spins and fij represents the coupling between them. 
This is invariant under rotations, yet in the ground state the spins are aligned 
giving rise to a magnetic field. So clearly in the ground state the symmetry 
of the Hamiltonian is spontaneously broken. 



Goldstone's theorem states that for a spontaneous breakdown of a symmetry 
there is an associated massless mode the Goldstone boson. In the case of 
QCD the axial symmetry is spontaneously broken. In this case Goldstone's 
theorem tells us that there should be a massless boson with the same quantum 
numbers as the as the generator of the broken symmetry. 




(4.25) 



4.4 Goldstone's Theorem 
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4.5 The relevant energy scale of QCD. 

In QCD the light quarks (up and down) have masses of approximately 5 
and 10 MeV. The relevant energy scale of the theory is given by Aqcd which 
is approximately 200 MeV. This is far greater than the light quark masses. 
Hence it would be reasonable to assume that the axial symmetry would be an 
approximate symmetry of QCD and so the axial current should be partially 
conserved. 



4.6 The PCAC relation 

Following on from the previous section we now investigate the current asso- 
ciated with the axial transformation in QCD. To do this we begin by consid- 
ering the weak decay of the pion. This is described by the matrix element of 
the axial current between the vacuum and pion states, i.e. 

(oKmAq)) (4.26) 

as before is the axial current and |7r^) is a pion state. By Lorentz symmetry 
this matrix element must be proportional to the pion momentum (g^): 

{0\A;m^\q)) = iUq,5^' 
where (O|0"(O)|7r^(g)) = 5"^ (4.27) 

/tt is a constant of proportionality (the pion decay constant) and is deter- 
mined experimentally. 5"'^ is used to normalise the states, and represents 
the pion fields. We take the divergence of the above equation and make use 
of integration by parts, we find: 

(0|a'^^^(0)|7r^(g)) = UmlS^'' 

= Uml{0\rm7r''{q)) (4.28) 
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This motivates the effective relation: 

= Umlr (4.29) 

In the hmit of vanishing pion mass the axial symmetry is exact. Since the 
pion mass is small compared to other hadronic masses we conclude the ax- 
ial current is approximately conserved and that the axial symmetry is an 
approximate symmetry of QCD. This is the PCAC relation. 



4.7 The Goldberger-Treiman relation 

We begin our discussion of the Goldberger-Treiman relation by considering 
the axial current for a nucleon: 

A'^^N = 9a^Nlt.l^—^N (4.30) 

Here ■j/'Ar is an isospinor representing the proton and neutron, = 1-25 is a 
constant that arises from the fact that the axial current is renormalised by 
25% as seen in the weak /3— decay of the neutron. Since the nucleon"^ has a 
relatively large mass we do not expect its axial current to be conserved. If 
we use the free-Dirac equation we can show: 

d^Al^^ = igMN^N^'r'^iJN + (4.31) 

Since the nucleon interacts strongly with the pion we assume that the total 
axial current is the sum of the pion and neutron contribution. 

Al = gai'Nl^lJ-^^N + Ud^r (4.32) 



■^We use the term nucleon to refer to either the proton or neutron. 
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Where we have used the PCAC relation (sec: 4.6). This current is only 
conserved if d'^A'J^ = so using equation 4.31, we find: 

d'^d.r = -gci^i'Ml^r'^'il^N (4.33) 

This is simply the Klein-Gordon equation for a massless boson coupled to a 
nucleon. Hence if we require the conservation of the total axial current (eq. 
4.32) for our system the pion must be massless. This is exactly what PCAC 
(sec. 4.6) predicts. By allowing for a pion mass we arrive at the modified 
Klein-Gordon equation: 

(a^9^ + rnDr = -gai^^Nl.r'^^N (4.34) 

By re- writing the coupling we have arrive at the Goldberger-Treiman relation: 

9.NN = 9a^ ~ 12.9 (4.35) 

Jtt 

The experimental value of the pion-nucleon coupling, which we will denote 
from here on, is g'^'^ ~ 13.2. This agreement is quite remarkable when we 
consider that we are relating the strong interaction of pions and nucleons to 
the pion decay constant and the nucleon renormalisation constant which are 
taken from the weak interactions. This relation will become important when 
we attempt to construct a chirally invariant Lagrangian. 



4.8 The spontaneous breakdown of the axial 
symmetry 

We appear to have found a contradiction, the mesonic spectrum does not 
appear to refiect the axial symmetry because of the mass differences between 
mesons (sec. 4.2). However, as we have seen (sec. 4.6) the weak pion decay 
indicates that the axial current is partially conserved (PCAC). As previously 
hinted at this problem can be resolved if the axial symmetry is spontaneously 
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broken. To see how this can occur we will consider the simple example of 
the Lagrangian (Ccsf) for a complex scalar field (0). 

Ccsf = \d,M'-Vm) (4.36) 
If we express the field in terms of its modulus (A) and phase {6) 




(4.37) 



it is easy to see that the potential (V^(|0|)) in the Lagrangian only depends 
on the modulus of the field. This can give rise to the classic Mexican hat or 
wine bottle potential. For such a potential, each point along the minimum 
has the same modulus but corresponds to a different value for the phase {6) . 
If we assume that the ground state selects the value ^ = for its phase, and 
if A = Aq we can expand about the ground state point. Using A = Aq + «, 
we find the Lagrangian reduces to: 

^csf = lid.af + lAlid.ey - y{^)-\^'y'i^)+- ■ ■ (4-38) 
Hence 

2 1 2t/// f 

m„ = -a V —= 

ml = (4.39) 

So we see rotations in the 6 direction correspond to massless excitations. 
The idea presented here will be useful when we study sigma models in the 
following section. 
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4.9 Sigma models and chiral invariance 
4.9.1 The linear sigma model 

In this section we will study a very simple Lagrangian involving pions and 
nucleons known as the linear sigma model. This was first introduced in 
1960 by Gell-Mann & Levy [9]. Previously we introduced the idea of chiral 
transformations of the quark fields that correspond to the pion (vr) (sec: 
4.2). We can also perform these transformations on the scalar meson (a). 
Altogether we find: 



Ay 
Aa 
Ay 
A^ 



a 
a 



cr 



(4.40) 



We see that although these may not be individually invariant under a chiral 
transformation the sum of their squares is: 



Hence, 



Ay 

Av 
Aa 



vr 



a 



TT 



vr 



2 - 2(1^7?'^ 



[TT + a ) ^ TT + cr 



(4.41) 



(4.42) 



We will use this to guide us while constructing our model. 



• We begin by considering the pion-nucleon interaction, this is described 
by a pseudo-scalar combination of the nucleon field multiplied by the 
pion field: 

gSi'l''r-^>- (4.43) 

where we have introduced the nucleon field ip. Under a chiral transfor- 
mation this transforms exactly as vr^ because the nucleon term has the 
same quantum numbers as the pion. For our condition for chiral in- 
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variance to be satisfied (eq. 4.42) we must have a term that transforms 
as a^. The simplest choice for this is: 

gMi^)(T (4.44) 

The sum of these terms (eqs. 4.43 & 4.44) gives the interaction term. 

We must account for the kinetic energy of the particles. For nucleons 
this is just the Lagrangian for free massless fermions and for the mesons 
we introduce an average for the a and tt— fields: 

Nucleons : iipc^ip 

Mesons: : ^{d^nd^n + d^nd^a) (4.45) 

We now need to introduce a nucleon mass term. Previously we showed 
that an explicit mass term breaks chiral invariance (sec: 4.1.2). The 
easiest way to introduce a nucleon mass term is via the coupling of 
the nucleon to the a— field (eq. 4.44). To do this we give the a— field 
a finite expectation value. Using the Goldberger-Treiman relation we 
find that the expectation value of the a— field has to be: 

{a) = U (4.46) 

This causes chiral symmetry to be spontaneously broken (see sec. 4.8). 
To give the a— field this expectation value we have to introduce a po- 
tential which has a minimum at a = 

Up until this point we have chosen terms that satisfy our chiral con- 
dition. But we have shown in section 4.1.2 that quark mass explicitly 
breaks chiral symmetry. So for our model to be consistent with na- 
ture we should include a small symmetry breaking term. We first write 
down a potential that is chirally invariant. A simple choice would be: 

V = V{n' + a') = ^{n' + a'-f:y (4.47) 
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To add a symmetry breaking term to this we recall that in QCD the 
axial symmetry is broken by a term which has the form SC^^^ = —niqq 
so it would be sensible to introduce a similar term, 5Csb = eo", with e 
being the symmetry breaking parameter. Doing this gives a modified 
potential: 

V = ^U^ + a'-u^y-ea (4.48) 



4 

Here we have introduced a general parameter uq. The only constraint 
that we place on this is that in the limit of e — then I'o ^ fn- 

Putting these terms together we find the Lagrangian for the linear sigma 
model: 



A / \ 2 

-^^(z^tV^tt" + ^ipa) - -[n'^ + a'^- ulj -ea (4.49) 

4.9.2 Properties of the linear sigma model 

We will now briefiy review the properties of our linear sigma model. 

• To preserve the Goldberger-Treiman relation our potential must have 
its minimum at for this to be the case the parameter to leading 
order in e is: 

^0 = - (4.50) 

• We can calculate the mass of the a particle by comparing our La- 
grangian with the Lagrangian for the Klein-Gordon equation, and re- 
calling {L = T-V). 

Ckg = i(9,$a^$) - im|^$^ 

-Q^ = <B (4.51) 



CHAPTER 4. AN INTRODUCTION TO CHIPT 



So for our Lagrangian we have: 



mt 



2A^^ + f (4.52) 



Our pion acquires a mass even though we have not exphcitly written 
one into our model as we did for the a particle. 



2 «9V 

m. — 



= f 7^ (4.53) 

(TO ■''^ 



We see that this fixes the symmetry breaking parameter (e) via e = fnml 
and hence we note that the pion mass-squared is directly proportional 
to the symmetry breaking parameter. 

Since o"o = fn the nucleon mass has a contribution from the explicit 
symmetry breaking. To see this we must split the nucleon mass into 
a contribution from the symmetric part of the potential and one from 
the symmetry breaking term: 



= g^cTo = i^uo + (4.54) 

The contribution from the symmetry breaking term is often referred to 
as the pion-nucleon sigma term (Stttv). Using equations 4.52 & 4.53 we 
may express this in terms of the pion and sigma masses: 

2 

6 TTl 
^nN = ^^TTVTi - gnfn^ (4.55) 

This term can be determined via the extrapolation of low energy pion- 
nucleon scattering data and is believed to be approximately 40 MeV. 

In section 4.9.1 we discussed adding a symmetry breaking parameter (e) 
to a chirally invariant Lagrangian. We required that this parameter had 
the same axial symmetry breaking properties as a mass term in QCD. 
With this in mind we can adjust its strength so that it reproduces the 
ground state properties of QCD, i.e. the pion mass. Hence it would be 
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reasonable to expect the vacuum expectation value of these symmetry 
breaking terms to be equal: 

{0\ea\0) = -(0|mgg|0) (4.56) 

If we recall that e = m^/^ and that (a) = (0|cr|0) = and also write 
out the average quark mass explicitly we arrive at the Gell-Mann- 
Oakes-Renner relation: 

ralf^ = -!Ihi±II^(^0\uu + dd\0) (4.57) 



4.9.3 The Non-linear sigma model 

The linear sigma model has one fundamental flaw, and that is that the 
(J— field cannot be identified with any physical particle. We can however 
remove the dynamical effects of this particle by sending its mass to infinity. 
This is done by assuming that the coupling (A) in our linear sigma model is 
infinitely large. The effect of this is to give the potential an infinitely steep 
gradient in the sigma direction. This causes the dynamics of our model to 
be confined to the minimum for the potential. This is a circle and is often 
referred to as the chiral circle it is defined by: 

+ vr^ = fl (4.58) 

We may now express the a and tt fields in terms of angles ($): 

(t(x) 

Where $ 

Hence to leading order the angles, $, can be identified as the pion fields. We 
immediately see that this satisfies the constraint of equation 4.58. The fields 



/ttCOS 



<i>(x) 



/.$sin(^ 



In + Oi<t>') 



v$*$ & $ = — 



(4.59) 
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may also be expressed using a complex notation: 



U (x) = exp yi 
= cos - 




) 



~17 



) 



+ zr"$ sin 





— (a + ir^n' 



.a 



(4.60) 



Here U{x) can be identified as a unitary (n/ x Uf) matrix (in our case U{x) 
is a 2 X 2 matrix) and is often referred to as the chiral field. Taking the trace 
of a combination of these matrices meets the constraint imposed by equation 



Because chiral symmetry corresponds to a symmetry with respect to a rota- 
tion around the chiral circle all structures of the following form are invariant: 



As we shall see this has non-trivial implications. 



4.10 The Weinberg Lagrangian 

To write down a Lagrangian for the non-linear sigma model in Weinberg's 
form we must use our findings from the previous section and redefine the 
nucleon fields: 



4.58: 





tr{U^U) trid^U^d^'U) 



(4.62) 



iIjw = AV' 



(4.63) 



With 




(4.64) 



This can be visualised as a dressing of the nucleon fields. We now have 
nucleon quasi-particles surrounded by a cloud of interacting mesons. 
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The pion-nucleon interaction term is now given by: 



cos( + r"*! sin ^-5- 



-gnfnip 



expl Z75- 



-MNtjjwiJw (4.65) 



Where in the last hne we have used the Goldberger-Treiman relation 
(eq. 4.35). Notice that using the redefined fields we have reduced the 
interaction term to the nucleon mass term. 

The kinetic energy term now becomes: 

Nucleons : iipc^ip = i%l}wl\^$h^ipw 

Mesons : \{^^,^x^^''K + ^^a^^'a) = ^tr{d^Wd^U) 

The A parameter has a spatial dependency through the $ fields and 
consequently the derivative in the nucleon term also acts on A giving 
rise to additional terms. We re-express this term as: 

^^/^vyAt^AtV^i^ = + ^f^V^ + r^,A^)tlJw (4.67) 

Where: 

and ^2 = U{x) (4.68) 

A^ and are the axial and vector quantities respectively. 

We need not consider the transformation of the chirally symmetric 
potential (eq. 4.47) because here dynamics are constrained to the chiral 
circle and on it the potential vanishes. 
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Using the above we can now write a complete Lagrangian for the non- 
hnear sigma model in the Weinberg form: 

Cw = i'w{i$ + YV^ + 7^75^M - MN)i'w + ^-^trid.U^d^'U) (4.69) 

It will prove instructive to expand the Lagrangian for small ^ / f-K ^ 1 fluc- 
tuations around the ground state. 

T^(^w/7^r«V'i4^)-($" X {d^^'^)) (4.70) 

J n 

4.11 Properties of the non- linear sigma model 

Here we briefly review the important properties of the Weinberg Lagrangian. 

• We clearly see the Lagrangian has a non-linear dependence on the 
$ -fields. 

• We have removed the unphysical a— field. 

• The coupling between the pions and nucleons now has a pseudo-vector 
form. It involves the derivatives of the tt— field, which we associate 
with the momentum of the tt— field. Along with this we now have an 
iso-vector coupling term. 

• We recall that at the expanded level the $— field can be identified with 
the TT— field (eq. 4.59). Hence explicit chiral symmetry breaking is 
introduced into our expanded Lagrangian (eq. 4.70) by an exphcit 
pion mass term. 

The final point to make in this section is that expanding the full La- 
grangian (eq. 4.69) to higher orders in the $— fields gives rise to extra terms 
which correspond to higher order corrections. These can be identified with 
trees, loops etc. 
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How to deal with these corrections in an ordered fashion will be the 
subject of our next section. 



4.12 Chiral Perturbation Theory 
4.12.1 The philosophy. . . 

The underpinning of chiral perturbation theory lies in a theorem of Wein- 
berg's. Generally speaking his theorem states that the most general effective 
Lagrangian will contain an infinite number of terms that satisfy the symme- 
try of the theory with an infinite number of free parameters. To make this 
a practical proposition we must have a scheme that tells us how to organise 
the terms and then assess the importance of the diagrams that are generated 
by the interaction terms from a given Lagrangian. This is the job of Chiral 
Perturbation Theory (ChiPT). The essential idea behind ChiPT is is to re- 
alise that at low energies the dynamics of the strong interaction should be 
dominated by the lightest particles of the theory (the pions) and the sym- 
metries of the theory (chiral symmetry). Hence physical processes should 
be expandable in terms of the pion's mass and momentum in a way that is 
consistent with chiral symmetry. Our goal therefore is to build a effective 
Lagrangian of the form 

Ceff = C2 + C4 + Cq + . . . 

00 

= 5^ An (4.71) 

n=l 

the subscripts refer to the order in momentum (i.e. the number of derivative 
terms) or the level of chiral symmetry breaking (i.e. £2 has one power of 
chiral symmetry breaking, m^). We note that to a given order the effective 
Lagrangian obtained from ChiPT should be consistent with QCD. We also 
note that ChiPT is not a perturbation theory in the usual sense, we do not 
expand in powers of a coupling constant. 
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4.12.2 Counting schemes 

To begin with we consider only the pion-pion interaction. A chirally invariant 
Lagrangian must be constructed using structures of the form: 



Also each chiral field {U) contains any power of the field, which can 
give rise to higher order diagrams. So to identify which structures we should 
include in our effective Lagrangian and by how much to expand each structure 
we count the powers of pion momentum that contribute to the process that 
we wish to study. 

4.12.3 Building a counting rule 

We will consider an arbitrary Feynman diagram that contributes to a scat- 
tering amplitude. 

• The diagram will contain a certain number of loops, which we will call 
L. 

• It will have a number of vertexes, call these V2n- 

• Each vertex will involve derivatives of the pion fields which we will call 



2n. 

• The number of internal lines associated with the vertex will be called 



We now need to determine the power (D) of the pion momentum (g) that 
the diagram will have (q^)- We must consider three points: 



tr{^^W^^'U) 

tr{^^W^^'U)tr{^^W^''U) 
tr[{df,Wd>'Uf] 



(4.72) 



1. Each loop involves a integration over loop momentum ~ q"^ 
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2. Each internal line corresponds to a pion propagator and hence carries 
momentum ~ X 

3. Each vertex involving derivatives of the pion field will contribute ~ g^" 

Using these points we can now write down the total power [D] of the 
momentum for the diagram under consideration: 

D = AL-2I + 2nV2n (4.73) 

71=1 

This can be further simplified by using a relationship between the number of 
Loops and internal lines and vertices of a diagram: 

^ = /-$^V2„ + l (4.74) 

Giving the simplified result: 

D = 2 + 2L + ^y2n(2n-2) (4.75) 

71 = 1 

We now have a scheme that tells us to which order of the expansion a given 
diagram will contribute. 

4.12.4 Obtaining an effective Lagrangian 

We now use the results from our previous section, where we considered the 
simple case of pion-pion scattering to write down an effective Lagrangian 
for this process. We begin by noting that the simplest chirally invariant 
combination of the chiral fields f/^f/ does not make any contribution to the 
dynamics because U'^U = 1. Hence the most simple contribution is given by: 

£2 = ^trid^U^d^'U) (4.76) 

note the subscript denotes the number of derivatives involved. Because we 
are considering pion-pion scattering we must expand to fourth order in the 
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pion fields: 

A = Ud,^)' + 7T^[($5^$)' - + (4.77) 

The first term can be identified with a free pion in the chiral limit and so 
it is the second term that describes the interaction and although it has two 
parts both contributions have the same number of derivatives and so in terms 
of our power counting scheme should be considered as one vertex, hence at 
lowest order we only have one diagram which is the simple pole diagram. 
Using the counting scheme developed in the last section (eq. 4.75) we can 
determine the order of the pole diagram. We note that the vertex function 
involves two derivatives of the pion field and so is equal to one for n = 1 
and is zero for all other n, and there are no loops. With this information we 
calculate that the chiral dimension of the pole diagram is -D = 2. 



4.12.5 Moving away from the chiral hmit 

Until this point we have been working in the chiral limit. For our theory 
to be a realistic one we must introduce chiral symmetry breaking into our 
Lagrangian. This is done by including terms of the form: 

tr{U + U^) (4.78) 

The simplest symmetry breaking term that we can include is: 

= 4 - im2$2 ^ 0($4) (4.79) 

We see that to leading order this can be identified with a pion mass term 
(we again note that the constant term makes no contribution to the dynam- 
ics). As previously seen with the chirally invariant terms we can include 
many different terms involving the above symmetry breaking term into our 
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Lagrangian, e.g. 

tr{W + U) 

tr{d,Wd^U)tr{W + U) 
tr{df,Wd^'U)tr{d^Wd''U)tr{W + U) 
tr[{d^Wd>'Uf]tr{W + U) 

Hence again a counting scheme is called for. This time we must not only 
consider derivative terms but also pion masses. To do this we simply modify 
our previous scheme (eq. 4.75) so that the parameter 2n now counts not 
only the derivatives of the pion fields at a given vertex but also the pion 
masses. The lowest order effective Lagrangian is now given by: 

C2 = ^trid^U^d^U) + l^L_JLtr{U^ + U) (4.81) 

Now the subscript tells us we are working at two derivative order or one 
power of chiral symmetry breaking (m^). Expanding to the lowest order in 
the pion fields reproduces the Lagrangian for a free pion: 

£2 = \{d,^? - + 0($') (4.82) 

For the case of pion-pion scattering we expand to fourth order in the pion 
fields. The lowest order effective Lagrangian (£2) fhat includes chiral sym- 
metry breaking reduces to 

(4.83) 

The adjustable parameters in the Lagrangian are the pion's mass and its 
decay constant. These should be fixed by experiment. This Lagrangian 
could then be used to calculate pion-pion scattering lengths [13]. 
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4.13 The Adelaide Method 

4.13.1 Introduction 

In the following two chapters (5 & 6) we employ and expand the Adelaide 
method [12, 24, 28, 29, 30, 32] for chiral extrapolations. The Adelaide fitting 
anzatz has been designed to take into account the non-analytic behaviour 
that arises due to the spontaneous breaking of chiral symmetry in QCD. 
This spontaneous breaking of chiral symmetry ensures that there is no simple 
extrapolation of hadron masses in terms of quark masses. This prompts us 
to use an effective field theory (Chiral Perturbation Theory) to guide our 
extrapolations. 

4.13.2 The Adelaide Anzatz 

In QED when placing an electron in the vacuum we must account for a 
cloud of virtual positrons that will surround the electron. This process oc- 
curs because the vacuum is a polarisable medium, electron-positron pairs are 
constantly being created and destroyed. This process is known as screening. 
A similar process occurs in QCD due to pion loops. This process means that 
hadron interactions cannot be treated as point-like in effective field theories 
that model QCD. The Adelaide anzatz introduces a new parameter to Chiral 
Perturbation Theory (A) that takes account of the finite size of a hadron and 
its surrounding pion cloud. 

Dimensionally regulated Chiral Perturbation Theory allows hadronic masses 
to be expressed as expansions in powers of the pion mass, known as Chiral 
expansions. These expansions prove to be very poorly convergent. This is 
because Chiral Perturbation Theory is effective up to about 4:7: f^^ ^ 1 [GeV] 
[11]. The pion mass terms in the chiral expansion arise from integrating to 
infinite momentum. Hence we are integrating past the cut off of the effective 
theory. This leads to a divergent series in m". 

Another way of viewing this is to note that Chiral Perturbation Theory is 
an effective theory and is valid only when the mass terms are not too large. 
When the pion mass becomes too large, their Compton wavelength decreases 
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so that the pions begin to probe the internal structure of quarks and gluons 
inside the hadrons. Since (dimensionally regularised) Chiral Perturbation 
theory assumes these hadronic fields to be fundamental, an increasing (and 
ultimately infinite) number of counter terms are needed to mop up for this 
discrepancy. 

The Adelaide anzatz states that the expansion is more naturally expressed 
in terms of the size of the extended source (A) divided by the pion mass {m-j^) . 
Hence 

If > A 

• The Compton wavelength (A = h/mc) is smaller than the extended 
source. 

• Pion loops are suppressed by A/m^. 

• Hadron masses vary slowly and smoothly with quark mass. 
But if m^r < A 

• The Compton wavelength is greater than the extended source. 

• This is equivalent to trapping a particle in a box. Causing multi-particle 
systems to arise. 

• Gives rise to rapid non-linear variations with pion mass. 

• The uncertainty in energy gives rise to pair production. 

=^ This causes particles to undergo self interactions. 
4.13.3 Interaction Lagrangians 

The self interactions that particles experience give rise to a self-energy. To 
understand where the equations describing self-energy come from we will 
consider the lowest-order effective ttN Lagrangian £^jv [10]. Using an ef- 
fective field theory allows us to remove some of the complications of QCD. 
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An effective Lagrangian that is consistent with the symmetries of QCD can 
(in the low energy regime) have all high momentum interactions integrated 
out, leaving behind nucleons and Goldstone bosons as the only degrees of 
freedom. 

C^N = ^{iIp-mN + y t'^Ts^^) ^ (4.84) 

Here \& is the nucleon doublet, Ip represents the covariant derivative, is 
the nucleon mass taken in the chiral limit and qa is the axial-vector coupling 
constant again taken in the chiral limit, is a Hermitian quantity known 
as the vielbein and it is given by 

= i[u''dfj,u — udfj,u'^] (4.85) 

With u representing the square root of the chiral fields (see eq. 4.60). To 
find the interaction term for this Lagrangian we must expand the chiral fields. 
We recall the complex notation for the chiral fields (eq. 4.60) since for that 
representation u{x) is simply given by 

u{x) = expl i (4.86) 



2f; 

On expanding u and and substituting into the vielbein we find 

u, = ^ + 0($3) (4.87) 

When this is inserted into the Lagrangian Ct,-n we find the interaction La- 
grangian 

^^nt = -^^^7^5r"5;.$"vl/ (4.88) 

In chapter 6 we use a SU (3) Lagrangian as our starting point. For full QCD 
this would be'^ 

= itT{Bu ■ VB) + 2Dii{BS^'{u^, B}) + 2Ftr(5S^[M^, B]) 

-l-symmetry breaking terms due to quark masses (4.89) 



'''Note the explicit symmetry breaking terms due to quark masses have been dropped. 
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Here B = B'^X"' with A'^ being the Gell-Mann matrices, u represents the 
velocity of the baryon in heavy baryon chiral perturbation theory, VB give 
the covariant derivatives of the baryon fields and 5''^ define the spin operators 
[15]. What should be noted is that for an SU{3) flavour symmetry two 
independent interaction terms appear. In equation 4.89 these interaction 
terms have the coefficients D and F. The full QCD Lagrangian that is the 
starting point for the derivation of the self energy equations of chapter 6 
can be found in [17]. From here we make the appropriate extensions to 
partially-quenched QCD. The full QCD Lagrangian of [17] also includes a 
clear outline of the symmetry breaking terms which arise from the quark 
masses. Although here I have only sketched a brief outline, a very good and 
detailed review of this can be found in [16] along with example calculations. 

It should be noted that similar effective Lagrangians can be derived for 
vector mesons. An example for the partially quenched case can be found in 
[18]. It is the interaction terms from the Lagrangian in [18] that subsequent 
calculations for the vector meson self energies in chapter 5 are made from. 

4.13.4 Self energy integrals 

As an example of how the self energy integrals are arrived at, we will consider 
the pion loop contribution to the nucleon self energy. 

For the sake of simplicity we will work in the chiral limit. The free 
propagator Sp{p) is modified by the self energy S(p) described by figure 4.1 

K 



P P-k P 

Figure 4.1: Pion loop contribution to the nucleon self energy. 



62 



CHAPTER 4. AN INTRODUCTION TO CHIPT 



The propagator is given by 



iSpip) 



I 



(4.90) 



f — niN + zO 



1+ 



Here O"*" represents the usual Feynman prescription for a relativistic Greens 
function. We now read off^ the Feynman rule for the vertex of an incoming 
pion with four momentum q and isospin index a noting that this is exactly 
the term that appears in our interaction Lagrangian 



Integrating over the loop momentum leaves an integral that describes the 
pion loop contribution to the nucleon self energy. 



The self energy integrals in chapters (5 & 6) are expressed in a three dimen- 
sional form where the time component has been integrated out. 

4.13.5 Regularisation and Renormalisation 

In order to perform a calculation within the frame work of a Quantum Field 
theory we must regularise and renormalise divergent loop integrals. Reg- 
ularisation and renormalisation is the two step process of the removal of 
infinite divergences. Regularisation describes the process of quantifying the 
asymptotically divergent components of loop integrals. Renormalisation is 
the subsequent removal of these divergences such that the results are rendered 
finite and any dependence on the regularisation procedure is removed. 

There are numerous schemes for renormalisation, in this section we will 
prove an equivalence between the minimal subtraction (MS) scheme and 

^For a complete list of Feynman terms see Appendix A of [19] or for Feynman rules see 
section 4.4 of [16]. 




(4.91) 




(4.92) 



(4.93) 
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the finite range regularisation (FRR) scheme that is central to the Adelaide 
method. 

FRR is a central reason for the Adelaide methods success. It allows us to 
replace the poorly convergent Chiral expansions that dimensional regulari- 
sation gives us with a highly convergent series [12]. 

We continue by considering the leading order non-analytic (LNA) be- 
haviour of the nucleon in the heavy baryon limit. In this limit the four- 
momentum is factored into a velocity dependent part and a residual momen- 
tum part. A projection operator is then used to split the baryon field into 
large and small fields. We begin by noting the formal chiral expansion for 
the nucleon mass 



= ao + a2ml + Eat (4.94) 

The coefficients ao and a2 come from the bare nucleon propagator and 
its leading quark mass dependence. In the heavy baryon limit the pion loop 
contribution to the self-energy at LNA, S^r, is given by (see eq. 6.7) 



TT T 

^NN XttJ-tt 

3_ 2 

327r/2^^ 

TT 



- A;^ + mi 



TT 



The integration has been done in the relative integral (/^). It's in- 
frared behaviour gives the leading non-analytic correction to the nucleon 
mass. Isolating the pole from the divergent part of the integral (Jtt) gives 

U = - [ dk{k' -ml) + - [ dk-^ (4.96) 
71 J 71 J k'' + m^ 

Here the second integral is simply equal to m^. In the minimal subtrac- 
tion scheme we absorb the infinite contributions from the first integral into 
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renormalised coefficients in our chiral expansion 



rriN = Co + C2ml + Xnml 
Co = ao + Xn- / k^dk 

TV . 



C2 = a2-x.- I dk (4.97) 



TT 



We now consider the finite range regularisation method. Here we use the 
nucleon's finite size and physical form factor to motivate the introduction 
of a regulator function, u{k), that vanishes sufficiently fast as A; ^ 00. The 
relative integral now becomes 

= - / dkjj^u\k) (4.98) 

We will set our regulator function to the simplest function that we can 
imagine, a sharp cutoff that is defined by some scale (A). Our regulator 
becomes 

u\k) = Q{A-k) (4.99) 

Our integral now has the upper bound of A rather that infinity. We may 
integrate this explicitly to find 

2A3 2A 2 2 1 / A \ 

In = ^ ml + -mitan-U — ) 4.100 

Sir TT vr Vm^/ 

Taylor expanding about the chiral limit gives 

2A3 2A , o 2 4 X 
U = - ml + ml-—ml + ... (4.101) 

Sit TT ttA 

We absorb these contributions from the integral into renormalised coefficients 
to give (eq. 4.102). 

2 

rriN = Co + C2ml + x^ml - X^ml— + ... 

^ _ 2A3 
Co — ao + Xtt-;: — 
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C2 = a2-Xn— (4.102) 

71 

This demonstrates a mathematical equivalence between the finite range reg- 
ularisation scheme and minimal subtraction since by sending the scale (A) 
to infinity we recover the minimal subtraction equations. Further discussion 
of this equivalence can be found in [20] & [12]. 
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Chapter 5 



An analysis of the vector meson 
spectrum 

5.1 Introduction 

In this chapter we apply the chiral extrapolation technique developed by 
the Adelaide group. It is designed to extrapolate lattice Monte-Carlo data 
using a finite range regulator prescription [12, 24, 29]. The following section 
lists the finite-range regulator forms for the self-energy of the p meson in 
the pseudo-quenched case. The derivation of this can be found in [32]. We 
use the Adelaide expressions for the self energy to fit data generated by the 
CP-PACS Collaboration [22] in section 5.3. Section 5.4 then gives details of 
the chiral fits. We then discuss varying the quantity used to set the lattice 
spacing in section 5.5.1. Finally we make predictions for the p, K* and </> 
masses along with predictions for the J quantity [21] and compare these with 
experimental results. 

5.2 The partially quenched ansatz 

In this section we study the form for the self energies and S^^ corre- 
sponding to Eqs. (3 & 4) in [24]. The processes responsible for these self 
energies are depicted in figure 5.1. 
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Figure 5.1: The first diagram gives rise to the leading non- analytic contribu- 
tion to the p self energy. The second diagram gives rise to the next-to-leading 
non-analytic contribution to the p self energy. The last diagram discribes the 
rj' contribution to the p self energy. These diagrams give rise to equations 
5.3, 5.2 & 5.4 respectivly. 
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Here though we consider the "pseudo-quenched" case, where valence and 
sea quarks are not necessarily degenerate. In [24] the case of full QCD was 
considered. We also consider the self energy contributions due to the double 
hairpin (DHP) diagrams. Our analysis is restricted to the case where the 
valence quarks in the vector meson are degenerate, i.e. k^^j 



K-2 



Throughout this chapter we will use the following notation. 
Mps(v){P, f^'sea] ^vai' '^vai) refers to the pseudoscalar (vector) meson mass where 
the first two arguments refer to the sea parameters and the last two refer to 
the valence quark masses. We will also use the following shorthand notation: 

Mnon-deg ^ M(/?, /C.^a; ^^sea, ^val) 
M'^^a = M(/3, /tsea; /tval, /tval) 
j^umt _ M (^(3 , Ksea'-i f^sea.^ f^sea) • 



where the superscript unit refers to the unitary data where 



val — "-val 



deg refers to the "degenerate" data where kI^i = k^^i and these are not 
necessarily equal to Usea] non-deg refers to the non-degenerate case where 

'^val 7^ ^val "^^^C OUC of thcSC is CqUal to Ksea- 

The total self energy is given by: 
where the individual terms are given by: 



^TTTT - -7-5- / —717^^777^^ TTTTT l^-^J 



(F + (M7/-''^^)2)(A;2 + (M^f )2) 
x((M;°"-''^^)2-(M^f)2) 
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lipgl r k^u\k) dk 

3vr2/.' Jo (P + {Mps^yy 

x{{M^fy-{Mp^f) (5.4) 

with ulik) = + (M^o-'ie3)2 

and AM^p = M^"-'^'^^ - mJ'^^ 

We note that (uJnik) + AM^^p) > for all quark masses and nontrivial mo- 
mentum considered in the lattice analysis. The constants in these equations 
are given by gi^pn = 16 [GeV~^], fp-^^ = 6.028. fip & fi^ are the (physical) p 
and 71 masses respectively. We take g2 = 0.75 which is the preferred value of 
[18] and U = 3/32 [GeV]. 

We use a standard dipole form factor, which takes the form 



= (A2 + p)2 
Ujru;{k) = u{k) 



Unnik) = u{k)u ^{^iJ?p/A - nl) 
The self-energy equations are discretised using: 



i^y i^z 

with k^^y^^ = — 

(^-^^ x,y,z 

We would like the finite range regulator to regulate the effective field 
theory when k^, ky, kz tend to infinity. Of course, once any one of the kx, 
ky, kz are greater than, say, lOA the contribution to the integral is negligible. 
Hence, we would like the highest momentum in each direction to be just over 
lOA. So we use the following to calculate the maximum and minimum for i, 
j, k above: 
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AO A a 
.lOA a 



2tx 



N, 



where [. . .] is the integer part. We study a range of values of A which are 
chosen based on the value of A^^ = 630 [MeV] used in [24]. The value of A 
is highly constrained by the lightest data point in the My versus Mpg plot, 
and since the data used in [24] includes a much lighter point than in this 
study, we use its value of A to guide our choice. 

Figure 5.2 shows the various self-energy contributions, S^^, S^^ and ^£)hp 
as a function of Mp°g~'^'^^ (see Eqs. 5.2, 5.3 & 5.4) for the representa- 
tive Ksca) = (2.10,0.1382) dataset (sec. 5.3) with our preferred value 
of A = 650 [MeV] 




Figure 5.2: The self-energy contributions (see Eqs. 5.2, 5.3 & 5.4) versus 



M 



non—deg 
PS 



data for the ensemble (/3, Kgca) = (2.10,0.1382). 



In section 5.4.3 we perform a highly constrained fit to the complete de- 
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generate dataset. We use this method to determine an estimate of the correct 
value of the A parameter. 

5.2.1 Double Hairpin Diagrams 

In eq. 5.4, T,^^^p, the double hairpin contribution to the vector meson masses 
self-energy is defined. This section sketches the derivation of this term [28, 



The rj' can occur as an intermediate state in the vector meson's self-energy. 
However, because of the special nature of the rj', care must be taken in order 
to correctly account for its propagator. In figure 5.3 the rj' intermediate state 
in the vector meson's self-energy diagram is represented. These diagrams are 
known as the "double hairpin (DHP) diagrams".^ The left-hand diagram 
has no sea-quark loop insertions in the t]' propagator, whereas the right-hand 
diagram has one such loop. In quenched QCD, only the left-hand diagram is 
present, and in full (unitary) QCD, both are present, together with diagrams 
including an arbitrary number of sea-quark loops. In the pseudo-quenched 
case, the same diagrams are allowed as in the full QCD case, except that 
the quarks in the loops have a mass, ml^^, which is not equal to the valence 
quark mass, ml^^. 

Concentrating on the pseudo-quenched case, the 77' contribution to the 
DHP can be written as 



where k is the momentum carried by the rj', and /Iq is the coupling of the 
quark bilinears to eachother via the gluons in the sea. Note that the first 
term in eq. 5.6 corresponds to the case of no sea quark loop insertion (i.e. 
the quenched case), the second term to one sea quark loop insertion, etc. 
Resumming the series in eq. 5.6 and sending yUo ^ 00 (since it corresponds 

^It turns out that the "single hairpin diagram" does not contribute to the vector meson's 
self- energy. 



32]. 




(5.6) 



5.3. OVERVIEW OF CP-PACS DATA 



73 



to a mass scale which is much larger than the pion mass) gives 



(5.7) 



There is a second diagram where the i]' enters which involves neither a 
single nor double hairpin. Its contribution is 



1 



Combining eqs. 5.7 & 5.8, and rearranging gives 



(5.9) 



which leads to the DHP self-energy term in eq. 5.4. Written in this way it 
is trivial to see that in the full QCD case, where Mpg = Mp°g~'^'^^ = Mpg, 
the DHP contribution gives zero as expected, i.e. it only gives a non-zero 
contribution for the pseudo-quenched case. 





Figure 5.3: The first of these quark flow diagrams shows the double hairpin 
and the second shows the double hairpin with one sea quark loop insertions. 
In each case the t]' is the meson propagating along the top of the diagram. 



5.3 Overview of CP-PACS Data 

In [22], the CP-PACS collaboration published meson spectrum data from dy- 
namical simulations for mean-field improved Wilson fermions with improved 
gluons at four different /5 values. For each different j3 value there are four 
different Kgea values giving 16 independent ensembles. We summarise the 
lattice parameters used in table 5.1. 
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In figure 5.4 we plot the unitary (i.e. k^^j = k^^^ = Ksca) pseudoscalar 
mass against the lattice spacing, a^o for the 16 ensembles in table 5.1 (Note 
that (Mp5**)^ is a direct measure of the sea quark mass as outhned in section 
4.9.2). Also included (for reference) are the mass values of the physical 
pseudo-scalar mesons vr, K, "77^" . Note the large range of both a and m'^c^ in 
the simulations, and that the lattice spacing, a, is primarily determined by 
the P value rather than the m^gg^ value. 

The physical volume of the lattice was held fixed at La ~ 2.5 fm for the 
P = 1.80, 1.95 and 2.10, but the P = 2.20 ensemble had a slightly smaller 
physical volume. A study of finite volume effects due to this is beyond the 
scope of this work, and we treat all 16 ensembles on an equal footing. The 
mass ratio Mps/My is related to the mass of the sea quarks used and varies 
from 0.55 to 0.8. The lattice spacing a varies from around 0.09 to 0.28 fm. 
In our study we consider the two cases where the scale is set using ro [21] 
and the string tension a. 

For each of the 16 ensembles we consider five Kvai values. Hence a global 
treatment of the data set yields a total of 80 {My^, Mpg) data points in the 
analysis. 

We generate 1000 bootstrap clusters for all Mps and My data using a 
Gaussian distribution whose central value and FWHM are the same as the 
central values and errors published in the table XXI of [22]. 

Our errors are totally uncorrelated throughout - i.e. each My (/5, Ksca] i^lai) 
bootstrap cluster is uncorrelated with the corresponding Mps{P, ^sea; «^vai' ^vai) 
bootstrap cluster. Also the M(/5, Usea] ^^iai' '^vai) data is uncorrelated with the 
M(/3', Kgga^; K^^i, K^^i) data, and, furthermore, M(/?, Ksee,; nl^i, nl^i) data is un- 
correlated with the M(/3, Ksea; «^vai' ^vai) data. 

Hence we expect the statistical errors in our final results to be overesti- 
mates of the true error because we have not benefited from any cancellation 
of statistical errors which should occur when combining correlated data. It is 
possible to estimate the increase in our errors due to the fact that we do not 
maintain correlations as follows. The ratio Mp^g^ /My™^ listed in table 5.1 
is obtained from our bootstrap data. Comparing this with the Mp^^^/My"-^* 
data in table XXI of [22] (which benefits from the cancellation of correla- 
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tions), we can see that ignoring correlations increases the errors very roughly 
by 20%. It is reasonable to expect a similar increase in errors for other 
quantities we study. 

The lattice spacings 0^0,0- are found from table XII of [22] using ro = 0.49 
fm and ^/o' = 440 MeV. As in the case of the meson mass data we generated 
1000 bootstrap clusters with a Gaussian distribution. 

The action used in [22] is tree-level, rather than non-perturbatively im- 
proved and thus is presumed to have some residual lattice systematics of 
0{a). We fit the data assuming both 0{a) and C(a^) effects in sections. 
5.4.2 & 5.4.3, and are thus able to obtain continuum predictions. 



> 

^ 0-6 



> 0.4 



0.2 



0.05 



0.1 



0.15 

a [fm] 





0.2 



X 


P = 


1.80 


A 


|3 = 


1.95 


□ 


P = 


2.10 





|3 = 


2.20 


< 


\ 




► 


K 




T 


7C 





0.25 



I ~^ I 
I -T^ I 



Figure 5.4: A plot showing the range of sea quark mass [Mp^g^Y and lat- 
tice spacing, a^o), covered by the CP-PACS data as displayed in Table 5.1. 
(Mp™*)^ is the pseudo-scalar meson mass squared at the unitary point, i.e. 
where /tvai = ^sca)- The experimental points for the n^K and "r/s" mesons 
are also shown for reference. 
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Table 5.1: The lattice parameters of the CP-PACS simulation used in this 
data analysis taken from[22]. The superscript unit refers to the unitary data, 
i.e. where kI^i = k^^i = Kgea- Note that the errors reported in this table are 
obtained with our bootstrap ensembles (see sec.5.4.1). 
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5.4 Fitting Analysis 

5.4.1 Summary of Analysis Techniques 

Our method is centred on converting all masses into physical units prior to 
performing any extrapolations. An alternative to this would be to extrapolate 
dimensionless masses (i.e. values in lattice units) as in [22]. We believe that 
our method has the following advantages: 

• We can combine the data from different ensembles and treat it in a 
global manner. If we left the masses in dimensionless units, we could 
not combine data from different ensembles due to differing lattice spac- 
ings. 

• Dimensionful mass predictions from lattice simulations are effectively 
mass ratios, and so we expect some of the systematic (and statistical) 
errors to cancel, e.g. M'^^"'^"' = M* x a"^ = M*/M^ x M^'^^ where 
VL is the quantity used to set the lattice spacing, a, the superscripts 
expt refer to the dimensionless lattice mass estimate and experimental 
value respectively. 

We consider two different methods for setting the scale. These are deter- 
mining the lattice spacing from the string tension (o") and from the Sommer 
scale (ro). We find one method is better than the other. This is outlined 
in section 5.4.3. Table 5.1 lists values for a^^ and a^. We also consider the 
effects of using other quantities to set the scale. 

We compare the Adelaide method with a naive polynomial fit. Our fitting 
functions take the following form, for the Adelaide fits 

^(M;^^)2-Stot = a, + a2{Mplf + a,{Mplf + a,{Mplf (5.10) 
where T^tot is from Eq.(5.1), and for the naive polynomial fit 

M'^'' = ao + a2{Mplf + a^{MplY + a^{Mplf (5.11) 
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We divide these fits into two further subcategories. The first category 
includes the above fits and is referred to as "cubic" since they include cubic 
terms in the chiral expansion of m'^ca ^ i^Ps)'^- '^^^ second category is 
formed from fits with the coefficient set to zero in equations 5.10 & 5.11. 
We call this category "quadratic" . 

We note that the dominant functional form of My with (Mp^Y is linear 
for example see figure 5.5. This fact is exploited in the above fitting functions. 



This is why the Adelaide fit uses y {My^y — T^tot on the left hand side 
rather than [My^Y — Y^TOT which would be an equally valid chiral expansion. 
It follows for the above argument that we can expect the a„ coefficients to be 
small for n > 4, and this is in fact what we find. In the following subsection 
we fit to equations 5.10 & 5.11 for the 16 ensembles in Table 5.1 separately. 
We then consider a holistic approach where we combine the data from all 16 
ensembles and perform a single global fit. 

5.4.2 Individual ensemble fits 

We first consider an individual analysis of the meson spectrum. This is done 
by treating the 16 ensembles listed in table 5.1 separately. We perform fits to 
the five [My^^ ^Ps) data points available from each ensemble. The fitting 
functions used are the Adelaide (eq. 5.10) and the naive (eq. 5.11) fitting 
functions. We restrict our attention to quadratic (ag = 0) chiral fits because 
there are only five data points available for each analysis. We use tq to set 
the scale and the A parameter for the Adelaide fits is set to A = 650 [MeV] 
which is our preferred value (see Sec. 5. 4. 3). 

The results for the coefficients ao,2,4 which are obtained by fitting My 
against Mps using both the naive (eq. 5.11) and Adelaide (eq. 5.10) fitting 
functions are listed in table 5.2. The fact that the coefficients are small and 
in most cases poorly determined supports our decision to fit to the quadratic, 
rather than the cubic chiral extrapolation form. 

Another important point to note is that there is a level of agreement 
between the naive and Adelaide ao,2 coefficients, although their variation 
with Ksea. tcuds to bc different. 
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Table 5.2: The coefficients obtained from fitting My data against Mpg using 
both the naive and Adelaide fits (i.e. eqs.(5.11 & 5.10)) for each of the 16 
ensembles listed in Table 5.1. As discussed in the text we restrict these fits 
to quadratic rather than cubic chiral functions (i.e. ag = 0). The scale was 
set from Tq. 
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We give a representative example of these fits in figure 5.5 using the 
ensemble {P,Hsca) = (2.10,0.1382). This ensemble's {a,m'^^^) coordinates 



are closest to the physical point 



a, m'i 



(0, rriu^d) for ensembles with 



La ~ 2.5 [fm] (see fig.5.4). In figures 5.6 & 5.7 we investigate the correlation 




M„ ' [GeV'] 



Figure 5.5: A plot of My^'^ versus Mps data for the ensemble (/?, Kgea) = 
(2.10,0.1382) together with the results of the quadratic Adelaide (Eq.5.10) 

and naive (Eq.5.11) fits. M^"* is defined as M^"'' = 



{Mi^^y - Etot for 
the Adelaide fit (i.e. the L.H.S. of Eq.5.10 - note E^ot is negative). 



of the (oq, 02) coefficients for both the Adelaide and naive fits. As expected, 
as Co increases, 02 decreases. Both methods show this trend to some extent. 
The figures also indicate that there might be a systematic variation of ao,2 
with arg. To investigate this further we plot ao and a2 against a^p (for both 
the linear and Adelaide fits) in figures 5.8 & 5.9. We use these figures to 
motivate a continuum extrapolation of the form 

^ ^cont I ■\r individual ^ /'K 10^ 

^0,2 — ^0,2 +^0,2 '^r-o {^■^'^) 



We list the results of these fits in table 5.3. 
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Figure 5.6: A scatter plot of a2 against ao for the Adelaide fit showing their 
mutual correlation. 
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Figure 5.7: A scatter plot of 02 against for the naive fit showing their 
mutual correlation. 
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Figure 5.8: A continuum extrapolation of the oq coefficient obtained from 
both the Adelaide and naive fits Eq.(5.12). 



The values of table 5.3 confirms a statistically significant 

0{a) effect in the Cq coefficient but is absent from the 02 coefficient. (Note 
we could have also performed a ffi which involves for example a 0{a^) term. 
We investigate these ffiting forms in more detail in the following section.) 

We now investigate the possibility of the lattice meson spectrum having 
a sea quark dependency. This is done by plotting the coefficients ao,2 — 
j^uidzvzduai agaiust (l/Mp^**)^. The results of this are shown in figures 
5.10 & 5.11. (Recall the superscript unit refers to the unitary data k^^j = 
K^^i = Ksca-) This is done because from the usual PCAC relation, (Mp^**)^ oc 
'^sea- We chose to plot [l/Mp^g^Y as the x— coordinate rather than (Mp™*)^ 
because this allows us to plot the quenched point at (l/Mp^**)^ = rather 
than at infinity. It should be pointed out that the physical point corresponds 
to (1/M]^^^*)2 ~ 50, consequently it is some way from our data. We subtract 
-^mdtviduai in the y— coordinate of figures 5.10 & 5.11 in a hope that we 
will be left with the residual m^„„ effects. This is done because, as we have 
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Figure 5.9: A continuum extrapolation of the 02 coefficient obtained from 
both the Adelaide and naive fits Eq.(5.12). 
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Table 5.3: The coefficients obtained from the continuum extrapolation of 
both the naive and Adelaide ao,2 values from Table 5.2 using eq.(5.12). 
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seen, variations in lattice spacing dominate those in m'^^^. 

Figures 5.10 & 5.11 show that there are httle if any significant mj^^ depen- 
dencies in ao,2- Moreover hnear fits to (1/Mp™*)^ produces a gradient which 
is almost zero within errors for the a2 term. We conclude this subsection by 
noting that we have not observed any evidence of unquenching effects in the 
the data. 
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Figure 5.10: A chiral extrapolation of the cq— Xq""'™*"'""''^ coefficient obtained 
from both the linear and Adelaide fits. Also plotted is the quenched data 
point, see sec. 5.4.4. The scale was taken from r^. 



5.4.3 Global fits 

In this section we treat the degenerate data from the 16 different ensembles 
as a whole data set. Doing this produces a data set containing 80 points 
(16 ensembles with five {My^,Mpg) values in each). Our hope is that this 
larger data set will constrain the fits allowing us to fit to more complicated 
functional forms and also produce a highly constrained set of fit parameters 
•^0,2,... • Figure 5.12 is a graphical representation of the 80 degenerate CP- 
PACS data points where we have set the scale from rg. 
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Figure 5.11: A chiral extrapolation of the 02— X™'^™*'^""' a coefficient obtained 
from both the linear and Adelaide fits. Also plotted is the quenched data 
point, see sec. 5.4.4. The scale was taken from vq. 



If we are to treat the data as a whole data set it is very important to 
model the lattice artefacts correctly. Table 5.2 along with the discussion in 
the previous section indicates a variation amongst the Oq values with lattice 
spacing, but the 02 coefficient is approximately constant with lattice spacing. 
Also recall that the 04 coefficient was undetermined. Hence we believe that 
allowing for variation in the coefficient due to the lattice spacing will be 
sufficient to correct any significant lattice artifacts. 

We use the above to motivate the following fitting functions. We define a 
modified version of the Adelaide and naive fitting function based on equations 
5.10 & 5.11. 



+a^{MplY + a^{Mplf (5.13) 
^a,{MplY + a,{Mplf (5.14) 
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Mp/ [GeV'] 



Figure 5.12: A plot of the degenerate CP-PACS data set. We have set the 
scale using ar-g. 
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As in the individual analysis (sec 5.4.2) we refer to the above fits as 
"cubic", since they include the term oc m^. As above we also perform fits 
with ttQ set to zero, referring to these as "quadratic" . 

We include corrections for 0{a) and 0{a^) lattice spacing effects in the 
fitting functions 5.13 & 5.14. This is because the lattice action used is tree- 
level improved, and so we expect it to contain O(a^) errors, but as shown in 
section 5.4.2 there is also some residual 0{a) errors. 

We have studied fitting functions that include 0{a, a^) terms in the a2 
(and even 04) coefficients to try to uncover lattice spacing effects in the 
higher order coefficients but we have found that these fits are unstable. This 
confirms the findings of our individual analysis reinforcing our belief that the 
discernible lattice spacing effects are contained in the Oq coefficient. 

In this section we study two different methods for setting the scale. We 
use both the Sommer scale (ro) and the string tension, (a). We summarise 
these different fits in table 5.4. In total we study 2^ fitting procedures, any 
one of these fitting procedures can be built by moving from left to right across 
table 5.4 and making a choice from the available options in each column. 



Approach 


Chiral Extrapolation 


Treatment of Lattice 


Lattice Spacing 






Spacing Artefact 's in ao 


set from 


Adelaide 


Cubic 


Qq term has 




i.e. eq.5.13 


i.e. 0{Mf,s) included 


0{a + a^) corrections 




Naive 


Quadratic 


ao term has 


a 


i.e. eq.5.14 


i.e. no C(M|,g) term 


only O(a^) corrections 





Table 5.4: The different fit types used in the global analysis. Fits for each of 
the 2^ choices depicted above were performed. 

We expect these fits to be highly constrained since they are performed 
using a data set containing 80 data points and the largest number of free 
parameters studied is six {qq, Xi, X2, 02, 04 and qq). 

When performing the Adelaide fits we must determine the correct value 
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of the A parameter (sec 5.2). This parameter is introduced in the Adelaide 
approach to model the size of the quasi-particle under consideration. It is 
this length scale that controls the chiral physics. Although it is not possible 
to allow A to be a free parameter in our fits we can derive the best value for 
A as follows. We manually vary the value of A and then plot the x^/d.o.f. 
as a function of A. Figure 5.13 is a graphical representation of this. When 
the scale is set from tq the value as a function of A exhibits the same 
functional form for all fits and they share a distinct minimum at A ^ 650 
[MeV]. Explicitly this means that the A parameter has no dependence on 
the order of the chiral expansion of our fits (i.e. expanding to 0{Mpg) or 
0{Mpg) has no effect on the correct value of A) also there is no dependence 
on our modelling of the lattice systematics in the oq coefficient (i.e. we can 
choose to use either 0{a + a?) or 0{a?)). 

When the scale is set using the string tension (a) figure 5.13 again shows 
that the all exhibit the same functional form. Now though we see that 
all fits share a distinct minimum at A ^ 550 [MeV]. The discrepancy in the 
correct value of A which arises when using different methods to set the scale 
will be addressed in section 5.5.1. However we see that now the 0{a + o?) 
fits give a far better that the 0{a^) fits. To investigate this further we 
have fitted the data using a fitting function with only 0{a) correction in the 

coefficient (i.e. eq. 5.13 with X2 = 0). The results showed near identical 

values as for the 0{a + a^) fits. This is indicative of a dominant 0{a) 
lattice-spacing systematic in the cases where the string tension is used to set 
the scale. We offer no explanation why this should be the case. 

These preferred values of A (550 [MeV] & 650 [MeV] for the a^- and a^p 
cases respectively) are used to perform the 16 global fits that are outlined in 
table 5.4. We list the results of these different fits along with the x^/d-o.f. 
in table 5.5. 

We now summarise the results of these of these fits (tb 5.5 and fig 5.13). 
• Fit approach 

The smallest x^/d-o.f. (and hence the best fit) is given by the Adelaide 
method. Moreover we see consistently smaller for a given fit using 
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Figure 5.13: A plot of /d.o.f against A. The dashed horizontal line repre- 
sents increasing from its minimum value by unity for the ro data (i.e. it 
represents one standard deviation), see sec. 5.5. The intercept of this dashed 
line with the curves (at A =630 and 690 MeV) is used to derive upper 
and lower bounds for the preferred A value. 



Fit Scale ag°"* Xi X2 as 04 ag X^/d.o.f. 
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Table 5.5: The results of the global fit analysis. Fits for all 2^ fit combinations depicted in table 5.4 are shown. 



5.4. FITTING ANALYSIS 



91 



the Adelaide method compared to the corresponding fit derived from 
our naive approach. This confirms that the Adelaide approach is the 
preferred chiral extrapolation procedure. 

• Chiral extrapolation 

In all cases the cubic chiral extrapolation (i.e. including a 0{Mpg) 
term) leads to a undetermined ag coefficient. We also observe the 04 
coefficient in the cubic fits becomes poorly determined compared to its 
quadratic chiral extrapolation counterpart. 

• Treatment of the lattice spacing systematics 

Studying the coefficients Oq and 02 in table 5.5, we see that in the case 
where rg is used to set the scale, the coefficients have little dependence 
on the type of lattice spacing correction used (i.e whether 0{a + a^) 
or O(a^) is used in the Oq coefficient). We do see a reduction in the 
error of the coefficient when only a O(a^) correction is used, this is 
most likely due to reducing the number of degrees of freedom. When 
the scale is set using a, we see that this is no longer true and that the 
coefficients Oq and 02 do depend on the treatment of the lattice spacing 
systematics. This supports the conjecture that setting the scale using 
the string tension leads to 0{a) systematics. 

• Setting the scale 

We see that in the Adelaide approach the is drastically reduced 
compared to when a is used to set the scale (fig 5.13). This along with 
the discussion regarding probable 0{a) systematics in the a data, give 
us reason to favour setting the scale using tq. In the case of the naive 
fits there is no clear preference between setting the scale from either ro 
or a. 



Using the above to guide our choice we select the quadratic chiral ex- 
trapolation method with C?(a^) corrections in the ao coefficient where the 
scale is set from ro to define the central value of both the Adelaide and naive 
fitting procedure. The spread from the other fitting types is used to define 
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the error. We make predictions for physical meson masses in section 5.5 for 
these fitting types. 

5.4.4 Quenched data. 

Along with the dynamical case, we have studied quenched data from [22]. 
This data was produced from simulations that use the same (gauge) lattice 
action as the dynamical case. We use values listed in table XIII of [22] for the 
string tension and tq to determine the lattice spacing where we take tq = 0.49 
[fm] and = 440 [MeV]. 

As in our dynamical analysis we represent the quenched data in [22] by a 
Gaussian distribution of 1000 bootstrap samples. We ensure the distribution 
has its mean equal to central value of the original data and its FWHM is 
equal to the error of the original data in [22]. Table 5.6 gives an overview of 
the parameters of these simulations. 

In [22] the quenched data was fitted with a linear fitting function of the 
following form (see eq.(59) and table XIV of [22]) 

aMv = A^ + B^{aMpsf. (5.15) 

This linear chiral fit is a simplified version of our naive fitting functions. It 
contains no 04 or ag coefficients (hence is linear in m^^"). To analyse this 
data we use the values of the coefficients Ay and (table XIV of [22]) of 
the fits performed in [22]. We do this because no individual masses for the 
quenched data are published. We convert these coefficients into dimensionful 
values via = /a and a2 = B^a. We study both cases where the scale is 
set using tq and the string tension. Table 5.6 lists the resulting dimensionful 
ao,2 coefficients. 

We investigate the dependency of these coefficients on lattice spacing by 
plotting the coefficients qq and 02 against lattice spacing a. Figures 5.14 and 
5.15 represent the case where the scale is set from tq. We see a clear lattice 
spacing dependency for the ao coefficient; the 02 coefficient also appears to 
exhibit a dependency on the lattice spacing. We model the lattice spacing 



2.187 
2.214 
2.247 
2.281 
2.334 
2.416 
2.456 
2.487 
2.528 
2.575 



Volume 



16^ 
16^ 
16^ 
163 
16^ 
243 
243 
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X 32 
X 32 
X 32 
X 32 
X 32 
X 48 
X 48 
X 48 
243 X 48 
243 X 48 



[fm] 



"0 



O2 



[GeV] [GeV- 



[fm] [GeV] [GeV-i] 



0.196t^ 

0.1871^ 

0.1751^ 

0.163±^ 

0.1490t?o 

0.12811i 

0.1201t^ 

0.1143^^ 

0.1072t^ 

0. 1003^3 



0.777^^2 
0.776ti^ 
0.787ti2 
0.824^1^ 

0.83511S 

0.860ti° 

0.8431« 

0.855t.i° 

0.857ti° 

0.859±?j 



0.421±1« 

0.42311? 

0.42611° 

0.391111 

0.3871^ 

0.37111J 

0.39411 

0.3841i« 

0.3851* 

0.38511 



0.208311^ 
0.198011^ 
0.18561?o 
0.17291^° 
0.1580+^ 



0.7331«o 

0.7321« 

0.7421« 

0.77811° 

0.787l«g 



0.4461° 

0.44811° 

0.45211? 

0.415111 

0.4111^ 



0.13611? O.8IOII; 0.39411? 



0.126811^ 
0.12081° 
0.11321?, 
0.1061? 



0.79811J 
0.81011J 
0.81211?, 
0.811? 



0.4161° 
0.40511J 

0.4061?o 
0.411^ 



:3 

Co 



Table 5.6: The lattice parameters of the quenched CP-PACS simulation used in this data analysis [22] together 
with the results of a linear chiral extrapolation. Note that the errors reported in this table are obtained with our 
bootstrap ensembles (see sec. 5. 4.1). 
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0.88 



0.86 



0.84 



0.82 - 




0.78 



0.76 



Figure 5.14: An O(a^) continuum extrapolation of the quenched coeffi- 
cients (i.e. using eq.(5.16) with x[^^ = 0). tq was used to set the scale. 



artefacts by assuming the same two a— dependencies in ao,2 as in sec. 5.4.3 

0(0,2) — 0(0,2) + ^1 a + y^2 o. [o.iQ) 

We perform the above 0{a + a^) continuum extrapolation, along with an 
(9(a^) extrapolation (i.e. we set xf'''^^ = in eq 5.16). We choose an 
0(0^) extrapolation rather than 0{a) because we expect the action to be 
dominated by O(a^) lattice spacing artifacts. Hence a linear extrapolation 
in the lattice spacing (achieved by setting X^''^^ = in eq 5.16) would not be 
appropriate^. Our results from these fits show that the 0{a + a^) fit leads to 
poorly determined coefficients xf'2^\ A point of inflection can also be seen 
in the fit which is located between the data and the continuum point (a = 0). 

We list our results for the values of a™;^* and X2°'^'' from the 0(0^) ex- 
trapolations in table 5.7. Graphically these fits are depicted in figures 5.14 

^Preliminary investigations into a linear continuum extrapolation lead to unreasonable 
continuum estimates for 00°^** 
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0.36 - 
O.35L 



0.01 



0.02 
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0.03 



0.04 



a [fm] 



Figure 5.15: An 0{a^) continuum extrapolation of the quenched 02 coeffi- 
cients (i.e. using eq.(5.16) with ' = 0). tq was used to set the scale. 



and 5.15. The fact that the X2 ' coefficients are significantly different from 
zero provide clear evidence for lattice spacing effects in the ao,2 coefficients 
for the quenched case. The continuum values, plotted in figs. 5.10 

& 5.11 and are surprisingly consistent with the dynamical data. 

In the following section we produce continuum estimates for the masses of 
the vector meson spectrum. We use these quenched values of aQ°2^ to include 
estimates of continuum quenched values of the vector meson spectrum 



5.5 Physical Predictions 

In this section we make physical predictions for the continuum masses of the 
p, K* k (j) {Mp.Mx* & M^). We do this for both the Adelaide and naive 
methods (eqs 5.10 & 5.11) that have previously been explored (sec 5.4.3). 

All mass predictions in this section are produced using our global fitting 
method rather than the individual analysis introduced in section 5.4.2. We 
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„cont. 

Uq 

[GeV] 


[GeV/fm] 


xl/d.o.f 


cont. 
2 

[GeV-i] 


[GeV^Vfm] 


xl/d.o.f 


From a^o 
From Oo- 


0.8951? 
0.85319 


-3.111 
-2.91^3 


6.7 / 8 
6.6 / 8 


0.3651? 
0.3811* 


1.51^3 
1.61^3 


9.9 / 8 
9.3 / 8 



Table 5.7: The coefficients obtained from performing a continuum extrapola- 
tion (eq. 5.16) to the ao,2 quenched coefficients. (Note we have set = 
- see text.) 



choose the global method because we expect the coefficients produced to be 
more accurate than those from section 5.4.2 since the global fits are highly 
constrained. 

In this section we also study the Mp prediction as a function of A. We 
choose to study the variation of the p-mass because it will be more sensitive 
to a changing value of A because S^ot is largest for smallest meson mass. 

The mass predictions for Mp, Mx* & M^ are obtained by setting Mp°^~'^^^ , 
Mpf & MPs'^ in equations 5.10 & 5.11 to the values outlined in table 5.8. 



Vector Meson M^f M^"""'^"^ M^f 



P 

K* 



f-^TT f^TT f^TV 

PK Pk/ Ptt 

Pr]s pK Ptt 



Table 5.8: Values for Mpf, Mpf'^"^ and M^f used in equations 5.1 - 5.4, 
5.10, 5.11 to calculate estimates of Mp, M^* and M,^. The symbol fi is used 
to represent the physical masses of the mesons. Note that the values for 
j^non-deg obtained by recalling that the non-degenerate meson contains 
one "valence" and one sea quark and that Mpg oc m^. 
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We assume an SU{3) flavour symmetry and so we expect the K and K* 
mesons to have valence quarks with a mass equal to half that of the strange 
quark. To calculate the self energy terms in the continuum a fourth order 
Runge-Kutta method is employed to calculate the integrals. We set Mp™* to 

throughout (because a unitary meson will comprise two sea quarks) and 
we assume that AM^p = 0. Physical predictions are made using all of the 2^ 
fitting types that were discussed in section 5.4.3 (table 5.4). To do this the 
coefficients, ao°"*&a2,4,65 from these fits (i.e. those in table 5.5), are used. We 
also make a prediction for the quenched vector meson spectrum as discussed 
in the previous section. Here we use the coefficients from section 5.4.4 listed 
in table 5.6. For all of the above cases we study both methods of setting the 
scale, using Tq and the string tension. 

We list results for all of our mass predictions in table 5.9. We have used 
our preferred values of A, A = 650(550) [MeV] (for the cases when the scale 
is set from ro(cr) respectively). 

Figure 5.16 graphically represents our investigations into how the Mp 
prediction varies with the value of A for each of the eight Adelaide fits. To 
estimate an acceptable range for the A parameter we use the plot of section 
5.4.3 (figure 5.13). Using this plot we can estimate the range of acceptable 
A values defined by increasing by unity from its minimum, this represents 
one standard deviation. 

The horizontal dashed line in figure 5.13 lies along values which are 
increased by one standard deviation (for our preferred method of setting the 
scale i.e the tq case). Hence we find that an acceptable range of values for A 
lies between 630 [MeV] < A < 690 [MeV]. 

We represent this range in figure 5.16 by plotting two vertical dashed 
lines at the acceptable maximum and minimum values of A. 

We summarise the information in table 5.9 and figure 5.16 below. 

• The statistical errors in the mass estimates are typically about 1%. 

• We see disagreement in the Adelaide fits when we choose to set the 
scale using different methods. The Adelaide procedure is very stable 
when we set the scale using tq. But when setting the scale using the 
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Quadratic chiral extrapolation contains 0{a + a^) 
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Quadratic chiral extrapolation contains 0{a'^) only 



Dynamical 
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0.77911 


0.8791^ 


1.020011« 


0.3891«g 


)7 


Adelaide 


a 


0.7741^3 


0.8531^ 


0.995011^ 


0.2991? 
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0.82511 
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0.4561« 
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Table 5.9: Estimates of Mp, M^^-., M,^ and J obtained from the global fits. 
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400 500 600 700 800 

A [MeV] 



Figure 5.16: A plot of Mp as a function of A from the Adelaide approach. 
Recall that the best A value when the scale is set from ro((T) is A = 650(550) 
MeV. The two vertical dashed lines define the range of acceptable A values 
(630 MeV < A < 690 MeV) obtained by increasing by unity in fig. 5.13. 
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string tension, we see that the four Adelaide fits do not agree and so 
appear to be unstable. We believe that this is most likely due to there 
being residual 0{a) errors when the scale is set using the string tension 
(sec 5.4.3). 

• The results obtained from the Adelaide fitting procedure are very ac- 
curate. At most they are twice the statistical standard error from the 
experimental value. For the Adelaide method to reproduce exactly the 
experimental Mp value, a re-adjustment of only around 1-2% in tq, and 
around 2-6% in ^/a^ would be required. 

• Notice (fig 5.16) that the variation of Mp with A is very small, it is 
about the same order as the other uncertainties. 

• The Adelaide method has central values that are far closer to the ex- 
perimental values than the naive method has. Furthermore the naive 
fitting method has larger spread of values than the Adelaide procedure. 

• The quenched results significantly overestimate Mp. The quenched 
value of the J-parameter is also significantly underestimated. (These 
two facts mean that the M^ quenched prediction is more accurate than 
the Mp value.) 

All of these points are in favour of the Adelaide method. Consequently 
we believe that the Adelaide method should be the favoured method when 
performing chiral extrapolations and we note that Adelaide method is a sig- 
nificant improvement over the naive approach. 

To give a final value for Mp for both the Adelaide method and the naive 
method, we use our preferred fitting function (the quadratic fit with O(a^) 
corrections in the ao coefficient) and our preferred method for setting the 
scale (from tq). Our error for the different fitting methods is obtained from 
the spread in the mass predictions (for the tq case). We also include an 
estimate of the error associated with the A parameter. We determine this by 
varying by unity (as described above). We can then simply read off this 
error from the vertical dashed lines in figure 5.16. 
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Hence Our final estimates are: 



M^'"^ = 779(4)tJ3+5jMeV] (5.17) 
^Naive ^ 825 (4)^12 [Me V] (5.18) 

(5.19) 

where the first error is statistical and the second is from the fit procedure. 
In the Adelaide case the third error is that associated with A. We do not 
make explicit any error that is associated with the determination of vq. 

We finally include a study of the J-parameter. This is normally defined 
as [21] 

(5.20) 



T dMv 
J = Mv 



dMpg K,K 

Here though we study the "discrete" version of this which we define as 

Tdiscrete ( — Mp 



J— = Mj,. [j^fzrj^ ) (5.21) 

We use this discrete version of J ( j'^^'^crete-j i^gc^^gg Iqq easily determined 
from experimental data, but J itself cannot. These two definitions coincide 
if My is a strict linear function of Mpg. 

Table 5.9 lists values for j'^^^crete^ ^^ig^t the value of J'^'*^^*^*^ is 

a severe underestimate of the experimental value. This is a well known 
phenomena and is no surprise. We also note that the estimates of J'^^'^^rete 
for the dynamical cases do increase toward the experimental value. 



5.5.1 Setting the lattice spacing 

In this section we investigate the differences in our results that occur when 
setting the scale using different methods. As previously mentioned we have 
studied two methods used for setting the scale. These are from the Sommer 
scale (tq) and from the string tension (cr). We investigate the ratio of these 
scales by plotting acr/arg against tq for each of the 16 ensembles in table 5.1. 
We see that this ratio is almost constant for the 16 ensembles and that there 
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is almost no evidence of 0{a) or rriq dependencies. The ratio is always greater 
than one and a rough estimate of its value would be around 5% above unity. 
We believe that this can be explained if the product v/a^o = 440 MeV x 0.49 
[fm] is approximately 5 % below its real value. It is this that presumably 
explains why, when setting the scale using tq, the estimates of the vector 
meson mass are larger than those found from using the string tension. Since 
the Adelaide method has a highly non-linear relationship between the lattice 
scale and estimates of the vector meson mass (My) due to the self 

energy (Etot) there is no corresponding simple relationship for the estimates 
of the vector meson mass made using different methods to set the scale. 
1.1 



1.08 



1.06 



1.04 



1.02 



0.98 - 



0.96 




Figure 5.17: A plot showing the ratio a„/arQ against a^o for the ensembles in 
table 5.1. 



This non-unit ratio must be responsible for the difference in the predicted 
best value of A that can be observed when using different methods to set the 
scale (fig. 5.17). It follows therefore if a^jar^ were unity, then we would see 
identical Adelaide predictions when using tq and a to set the scale. 

We also investigate one final method for setting the scale. This is the 
method of [25] . This method fixes the lattice spacing (and the strange quark 
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mass) from the {K, K*) mass point, i.e. using the J-parameter. The results 
of this method are graphically represented in figure 5.18. As in our global 
analysis we plot the full degenerate CP-PACS data set (80 degenerate points 
from [22]). Figure 5.18 also includes the unitary UKQCD points from [26]. 
It is quite remarkable that when this method is used, our data lies on an 
almost universal straight line. Compare this with the case where ro is used 
to set the scale (fig. 5.12) and note that this is exactly the same data set. 
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Figure 5.18: The vector meson masses. My versus Mp^ from the CP-PACS 
collaboration [22]. The scale is set from the (fC, K*) point using the method 
in [25]. Also shown are the unitary UKQCD points from [26]. 



The method of [25] seems to be an ideal way of setting the scale since the 
data would be well modelled by a simple linear fit. This is because the data 
is forced to go through the {Mk, Mk*) point, leaving the gradient as the 
only free parameter. But in doing this we are normalising away the expected 
non-linear behaviour as the chiral limit is approached, and it is exactly this 
behaviour that we try to describe using the self-energy term Stot in eq. . It 
is for this reason that we believe, despite the universal behaviour of the data 
when the scale is set using the J-parameter (fig. 5.12), it would be incorrect 
to use this method to set the scale in this work. 
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5.6 Conclusions 

We conclude this chapter by hsting the results of our study. 

• We have shown that the Adelaide method is a valid chiral extrapola- 
tion procedure and we have generalised the Adelaide chiral ansatz to 
"pseudo- quenched" case (i.e. when Kgea 7^ ^^vai)- 

• We have quantified the residual 0{a) effects in the CP-PACS data. 
(See e.g. figs.5.8 & 5.9.) 

• We have studied different fitting methods and found our global proce- 
dure to be the better method. 

• We have demonstrated that the Adelaide method can predict a pre- 
ferred value for the A parameter. 

• We have indicated that small errors in the values of and a might 
be the cause of the slight inaccuracy in the central values of our mass 
estimates. 

• We have obtained estimates of the p, K* & masses with tiny (statis- 
tical) error bars. (See table 5.9.) 

• We have estimated systematic errors in the p mass from the fitting 
procedure (both chiral and continuum fitting procedure). 

Note that we have not modelled finite-size effects - see Sec. 5.3. This is 
because we do not have enough different volumes to undertake such a study. 
Finite volume effects are considered by the Adelaide method (eq 5.5) since 
the momentum integral is replaced by the appropriate kinetic sum. 



Chapter 6 



An analysis of the Nucleon 
mass from lattice QCD 



6.1 Introduction 

In this chapter we again use the chiral extrapolation technique developed by 
the Adelaide group to estimate the mass of the nucleon from lattice QCD. 
We will employ a dipole form factor as we did in chapter 5 and we will 
also study a Gaussian form factor in an attempt to prove that the Adelaide 
method is not dependent on the finite-range regulator that is employed. In 
the next section we list the finite-range regulator form for the self-energy of 
the nucleon in the pseudo-quenched case. We again use the data generated by 
the CP-PACS group in [22] . We provided a comprehensive review of this data 
in 5.3 for the p case and also, for the nucleon data, in 6.3. In section 6.4 we 
outline the various fitting methods that we employ. Section 6.4.4 investigates 
the differences between the Gaussian and dipole finite-range regulators. The 
section following this contains our physical predictions for the nucleon mass. 
We then discuss the different methods of setting the scale. Finally in section 
6.7 we draw our conclusions. 
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6.2 The partially quenched ansatz 



In this analysis we restrict our attention to the case where the valence quarks 
are degenerate [kI^^ = k^^^ = K^^y = K^ai)- 

We begin by defining the following shorthand notation 

^non_aeg ^ ^^^^^ 

) sea 1 ^sea ; val ; val j 
; '^sca; ^vah ^vah ^valj 
j^non_deg ^ Mpcj /tgea; ^^sea, ^^val) (6-1) 

; '^sea; ^vah ^valj 
; '^^seai ^sea; '^^seaj 

where Mb{ps) is the Baryon (pseudo-scalar) mass with B = N A. The first 
two arguments of Mb{ps) refer to the sea structure (i.e. the gauge coupling 
and sea quark hopping parameter) and the last three(two) arguments refer 
to the valence quark hopping parameters. 

We also define the following integrals 

""^ k^u\k)dk 



I{Mps,6M) 



Tx Jq uj{uj + 6M) 



h{Mps) = - / ^ (6.2) 



Where we have used: 



u{k) = sJk^ + Mls (6.3) 

Here Mps can be Mp^ or Mp°^~'^'^^ . We define this along with values for 
6M explicitly in the individual self energy terms below. 

We use a standard dipole form factor, which takes the form 
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We also study a Gaussian form factor 

u{k) = exp(^-^^ (6.5) 

The self energy (Sat) is the total contribution from those pion loops which 
give rise to the LNA and NLNA terms in the self energy of the baryon, and 
also the contributions that arise from the rj' diagrams. Explicitly we write 
the processes as ^ iVvr ^ iV, iV ^ Avr ^ iV, iV ^ Nr]' N and 
N ^ At]' ^ N (figure 6.1). 



71 






Figure 6.1: The four diagrams that give rise to the leading and 
next-to-leading non-analytic contributions to the nucleon mass along with 
the DHP contributions from the f]'. This diagrams give rise to equations 6.7. 



In the limit of full QCD these rj' contributions vanish. For partially 
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QQCD in the heavy baryon hmit this may be expressed as [28]: 

= (Tnn + ^NN + ^NA + ^TVA (6-6) 

Exphcitly we have: 

_ 3(F + Df (j.^^deg Q^ 

^NN - 327r/2 \/^^^^PS^^) 

+a(/(M^^"-'^"^ ^non-des _ ^deg^ _ J (^M%' , 0)) 



'NN 



327r/2 



1 1 



3271 3' \ ^ ' ^ 

-M^"^) - /(M^f , M^*^^ - M^^^)) (6.7) 



The parameters a, /3 & 7 are derived from the standard SU{6) couphngs^ 
[29] exphcitly we take 

loops 

a -- 



P 



2{F + Df 
loops 



2{3F - L>)2 
7 = -2D 

loops = -{3F + Df + 3{D ~ F)^ (6.8) 
We use the constants F = 0.51 and D = 0.76 which are determined from 



-^For a full discussion see [31]. 
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fitting semi-leptonic decays at tree level e.g. [14]. 
Our fitting function takes the following form 

Mf;' = a, + a2{Mplf + a,{Mplf + a,{Mplf + T.^ (6.9) 

This equation is based on the chiral expansion in [29] and previous work in 
[30]. It also enjoyed considerable success in [32]. 

The self energy integrals are discretized using the same method outlined 
in section 5.2 of chapter 5. 

A value for the A parameter is determined by varying A and looking for 
a mimima in the /d.o.f . 

Figure 6.2 graphically represents the different contributions to the nucleon 
self-energy, along with the physical continuum values for pion processes. We 
note that there are no continuum values for the rj' since these processes 
disappear for physical values of the parameters (as required). 

6.3 The CP-PACS Nucleon data 

In this chapter we again use data published in [22]. Here we use the baryon 
data though. As before this data comes from dynamical simulations for 
mean-field improved Wilson fermions with improved gluons. We again study 
four different (3 values which each have four different Ksca values, giving 16 
independent ensembles. The lattice parameters used have been summarised 
in table 5.1. Figure 5.4 is a graphical representation of the unitary pseudo- 
scalar masses plotted against the lattice spacing a^o and we again recall that 
(Mp5**)^ is a direct measure of the sea quark mass (sec. 4.9.2). In this chapter 
we consider the same two methods of setting the scale as in chapter 5, namely 
from the string tension (a), and from the Sommer scale (ro). The degenerate 
data set contains 80 data points (five Kvai values for each (/3, Kgca) point), 
and as before we generate a 1000 bootstrap clusters for all Mps and 
data. The data has a Gaussian distribution with a central value equal the 
Mps{Mn) value published in table XXI(XXII - for the degenerate nucleon. 
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Figure 6.2: Here we plot the self-energy contributions (Eqs. 6.7) versus 
{Mp°^~'^^^Y for the entire degenerate data set (dashed lines are a guide for 
the eye only). We use the dipole form factor and choose an arbitrary value for 
the Lambda parameter, A = 1 [GeV]. We also include continuum data (the 
straight lines) for the pion processes (the eta case vanishes in the physical 
limit) which is obtained by solving the self energy equations using a fourth 
order Runge-Kutta algorithm. 
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XXIII - for the non-degenerate nucleon^) of [22] and has a FWHM equal to 
the published error. As before (sec. 5.3) we use totally uncorrelated data 
throughout, hence we expect our statistical errors to be overestimates of the 
true error in our results. As before the lattice spacings and a^g are taken 
from table XII of [22]. We generate 1000 bootstrap clusters with the correct 
Gaussian and FWHM distribution. The values tq = 0.49 and ^/o' = 440 MeV 
are used. Again we assume that the data has both 0{a) and O(a^) lattice 
systematics which we investigate in section 6.4. 



6.4 Fitting analysis 

6.4.1 Summary of analysis techniques 

The philosophy behind our fitting method remains the same as for our inves- 
tigation of the meson spectrum (chap. 5), i.e. we work in physical units when 
performing our extrapolations. We do this because it allows us to combine 
data from different ensembles; this cannot be done for the dimensionless data 
because of differing lattice spacings. Also we expect that we will benefit from 
some cancellation of the systematic (and statistical) errors. This is because 
dimensionful mass predictions from lattice simulations are effectively mass 
ratios (sec. 5.4.1). 

As mentioned previously we study two methods for setting the scale, using 
the string tension (a) and the Sommer scale (ro). As with the meson analysis 
(chap. 5) we find a preferred method for setting the scale which we discuss 
in section 6.4.3. 

We compare the Adelaide fitting procedure to the corresponding naive 
fitting function. After trying many different ways of fitting the data and 
allowing the error in the coefficients and the for these fits to guide us, we 
find that the data is best fitted by the following 



^We take the values of ms in table XXIII to be the mass values for the non-degenerate 
nucleon. We can do this since the interpolation operator for N and E have the same 
quantum numbers 
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Mn-^n = a^ + a2{Mplf + a^{Mplf + a^{Mplf (6.10) 
The corresponding naive fitting function is 

= ao + a2{Mplf + ai{Mplf + a^{Mplf (6.11) 

Again we divide these fits into two categories. These are referred to as 
"cubic" and "quadratic" . We remind the reader that this is because "cubic" 
fits include cubic terms in the chiral expansion of ml^^ oc (Mp*^)^. The 
"quadratic" fits have the ae term set to zero and so have terms that are 
quadratic in ml^^. We choose to use a subtracted fit for convenience. A 
fit that has the self-energy added to the RHS of equation 6.11 would be 
equally valid. In figure 6.3 we plot the dimensionful nucleon data and also 
the subtracted nucleon data. We see that the original data for the nucleon 
has some curvature but this is mostly corrected by subtracting the self-energy 
from it. So for high values of A we expect the higher order coefficients in the 
cubic fits to be poorly determined or approximately zero for the Adelaide 
case. Although this is the case, we will see that the data is better fitted by 
a moderate value of A. 

In the next subsection we fit to equations 6.10 and 6.11 for the sixteen 
individual ensembles (table 5.1). Following this we fit to 6.10 and 6.11 for 
the entire degenerate data set (as in sec 5.4.3). 

6.4.2 Individual ensemble fits 

In this section we treat the sixteen ensembles separately. We do this by fitting 
to the five degenerate data points {M'^^ , Mp^) in each (3, Kgoa ensemble. We 
use our preferred form factor which is the dipole form factor (sec 6.4.3 & sec 
6.4.4) and we only consider the case where tq is used to set the scale. For 
the dipole form factor we choose a preferred value for A (when setting the 
scale from tq) of 600 [MeV]. Section 6.4.3 shows that the preferred value of A 
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Figure 6.3: The nucleoli mass and the subtracted nucleon mass versus 
[MpgY for the entire degenerate data set. Again we use the dipole form 
factor and we use the same arbitrary value for the Lambda parameter, A = 1 
[GeV]. The scale is set using the string tension (a). 
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has a shght dependency on the modelhng of the lattice systematics, but this 
is a very small variation of approximately 1%. The fits considered in this 
section are quadratic, i.e. we set = 0. This is done because cubic fits for 
the individual case all have 100% error in the and ag coefficients. Table 
6.1 lists the coefficients for both the Adelaide fits and also the naive fits. 
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Table 6.1: The coefficients obtained from fitting Mn data against Mpg. We 
list results for both the naive and Adelaide fits (eqs. 6.11 & 6.10 respectively) 
for each of the 16 ensembles listed in Table 5.1. A dipole form factor was 
employed for the Adelaide fits using A = 600 [MeV] and the scale was set 
from tq. 

As expected the leading Adelaide coefficient is always greater than the 
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corresponding coefficient from the naive fits (oq'^'^' > Oq"™*^). In nearly all 
cases the 02 coefficient is smaller for the Adelaide ffis [a'^'^'^^ < a'^"'™'^). The 04 
coefficients are approximately the same for both ffis (a^'^'^' ~ aj"™*^), but the 
error in this coefficient is very large, typically 50%. We see only a few cases 
where the 04 coefficient is zero within errors though, indicating its presence 
is needed. Importantly nearly all of these coefficients are equal within errors. 

There appears to be no overall trend with the sea quark mass for any of 
the coefficients. This indicates that unquenching effects are minimal. 

Figure 6.4 is representative of all ffis. It comes from the Kgea) = 
(2.10,0.1382) ensemble. This data set is one of the closest to the physical 
point (fig 5.4). 




Mp/ [Gev'] 

Figure 6.4: A plot of Mj^ versus Mp^ for the ensemble (/9, Kgea) = 
(2.10,0.1382). Included are the results of the quadratic naive (Eq. 6.11) 
and the quadratic Adelaide (Eq. 6.10) ffis. The scale is set from tq, we use 
a dipole form factor and our preferred value for A (A = 600 [MeV]). 

We now go on to investigate the correlation between the coefficients Oq 
and a2. As can be observed in figures 6.5 and 6.6 there is a clear and well 
defined correlation between the qq and 02 coefficients. As expected, when 



116 



CHAPTER 6. AN ANALYSIS OF THE NUCLEON MASS 



the value of increases, the value of 02 decreases. The correlation here is 
far better than that observed in section 5.4.2 for the meson data. 
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Figure 6.5: A scatter plot of 02 against qq for the Adelaide fit investigating 
their mutual correlation. 



We now investigate the variation of the 00,2 coefficients with the lattice 
spacing 0^3. As in section 5.4.2 we plot 09,2 against a^p and use these plots 
(figs. 6.7 & 6.8) to motivate the following continuum extrapolation 

^ ^cont I \rindividual ^ ! R ^ n\ 

«0,2 — 00^2 +^0,2 ^r-o lO-J-^j 

Interestingly we investigated a continuum extrapolation of the form 
O'O = <^o°^^ + Xq'^'^™*'^""' but found that the data was better fitted^ by equa- 
tion 6.12. This is the case for both the Adelaide and naive data. It is at odds 
with section 6.4.3. The results of the fits corresponding to equation 6.12 are 
listed in table 6.2. 

We see that although errors are high for all values of X*"'^™^'^""' the better 

■^A reduction in the of about 5% was observed for the naive fit, and of around 6% 
for the Adelaide fit. 
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Figure 6.6: A scatter plot of a2 against ao for the naive fit investigating their 
mutual correlation. 
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Figure 6.7: A continuum extrapolation of the ao coefficient obtained from 
both the Adelaide and naive fits Eq.(6.12). 
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Figure 6.8: A continuum extrapolation of the 02 coefficient obtained from 
both the Adelaide and naive fits Eq.(6.12). 
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Table 6.2: The coefficients obtained from the continuum extrapolation of 
both the naive and Adelaide 09,2 values from Table 6.1 using eq. 6.12. 
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determined lattice spacing effect appears to be present in the leading coef- 
ficient (ao). Because errors are so high in these coefficients we believe that 
this is a sign of there being minimal 0{a) effects in the lattice data when the 
scale is set using a^o. We investigate this further in section 6.4.3. 

To close this section we provide a brief investigation of the lattice nucleon 
data having a sea quark dependency. To do this we use the same method 
as in 5.4.2 i.e. we plot ao,2 - X'J'f"''^''^^ against {M^f^f. We remind the 



reader that we subtract X, 



individual 
0,2 



ttro in the y— coordinate in an attempt to 



leave the residual m'^ei^ effects. The results of this can be seen in figures 6.9 
& 6.10. We observe no discernible trend with the sea quark mass (from the 
PCAC relation, (Mp^**)^ oc m'^^.g) for either of the coefficients. It can be 
seen however that as the mass of the pseudo-scalar decreases, the difference 
between the Adelaide and naive data points increases. 



1.25 



Adelaide data 
Adelaide fit 
naive data 
naive fit 




, unit^2 , ,-2-, 

(Mp„ ) [GeV ] 



Figure 6.9: A chiral extrapolation of the _ Xq^'^^^'^^^'o coefficient obtained 
from both the linear and Adelaide fits. The scale was taken from Tq. 



We conclude by noting that we have seen no strong evidence for un- 
quenching effects in the data. 
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Figure 6.10: A chiral extrapolation of the 02— Xg^'^^^'^^'^'a coefficient obtained 
from both the hnear and Adelaide fits. The scale was taken from tq. 
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6.4.3 Global fits 

We now analyse the complete degenerate data set by treating the 16 ensem- 
bles of section 6.4.2 in a global manner. This gives us 80 data points to 
work with. The larger data set should produce highly constrained fits (as 
seen in section 5.4.3). It is our hope that this larger data set will allow us to 
determine the higher order coefficients in our fitting functions (Eqs. 6.10 & 
6.11). We give a graphical representation of the degenerate (80 data points) 
CP-PACS nucleon data in figures 6.11 and 6.12 where the scale has been set 
using the Sommer scale (ro) and the string tension (cr) respectively. As can 
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I 

ll , \ , \ , \ , \ , \ 
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Figure 6.11: A plot of the degenerate CP-PACS nucleon data set. Here the 
scale is set using tq. 

be seen the data in these plots has clear curvature, moreover the agreement 
between data from different (/?, Usea) values is quite remarkable when we con- 
sider that these plots are representative of the degenerate data set and have 
had nothing done to them to correct for lattice artefacts. When treating the 
data in a global manner it is very important to ensure that the lattice spac- 
ing artefacts are modelled correctly. To this end we use the investigations 
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Figure 6.12: A plot of the degenerate CP-PACS nucleon data set. Here the 
scale is set using a. 



of the previous section and all we have learnt from the mesonic data (chap. 
5) to guide us when trying to account for lattice spacing effects. We have 
studied many different fitting functions where we try fitting forms that in- 
clude 0{a, a^) terms in the 02 and higher coefficients. We have found though 
that these fits are unstable. Hence we conclude that for the global method 
lattice spacing artifacts are dominant in the leading coefficient (oq). Hence 
we believe that accounting for 0{a) and O(a^) errors in the ao coefficient 
will be enough to correct for any significant lattice spacing effects present in 
the data. We choose 0{a) and O(a^) corrections because the lattice action is 
tree- level improved and so should contain O(a^) errors together with a small 
amount of 0{a) errors but as seen in chapter 5, when the scale is set using 
the string tension, 0{a) errors seem to dominate. Therefore our preferred 
fits, which are modified versions of equations 6.10 and 6.11, are 

Mn-Ej, = {ao + X^d') + a2{Mplf + a^{Mplf + a^{Mplf (6.13) 
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and the corresponding naive fitting function for the global case is 

= (ao + X„a")+a2(MS)' + 04(MS)' + a6(MS)' (6.14) 

As discussed in the summery of our fitting analysis (sec 6.4.1) we divide these 
fits into two categories. These are referred to as "cubic" and "quadratic" . As 
before this is because "cubic" fits include cubic terms in the chiral expansion 
of m^gjj oc [MpgY. The "quadratic" fits have the term set to zero and 
so have terms that are at most quadratic in ml^^. These two categories are 
divided into two further sub-categories. The first sub-category contains the 
fitting functions that have corrections for 0{a) lattice spacing effects in the 

coefficient, i.e. we have n = 1 in equations 6.13 and 6.14. The second 
sub-category contains the fitting functions that have O{o?) corrections for 
lattice spacing effects in the coefficient, i.e. we have n = 2 in equations 
6.13 and 6.14. Hence for the global fit analysis the maximum number of fit 
parameters in any one fitting function is five. Our data set contains 80 points 
so we hope that this method will provide highly constrained fit parameters 
compared to those from the individual fitting method (sec 6.4.2). This is 
because for the individual fitting method the number of fit parameters was 
three, but the number of data points in each (/3, Ksca) ensemble was only five. 

In this subsection we study two different methods for setting the scale. 
We remind the reader the previous section only studied one method of setting 
the scale. Here we set the scale using the Sommer scale (ro), and from the 
string tension, (a). Finally for the Adelaide method we have two different 
form factors to study (as outlined in sections 6.1 & 6.2). These are the 
dipole and Gaussian form factors given by equations 6.4 and 6.5 respectively. 
Hence for the naive fitting method we have 2^ different fitting procedures, 
and for the Adelaide method we have 2^ different fitting methods. These 
fitting methods are summarised in table 6.3. As in section 5.4.3, any one of 
these fits can be built by moving from left to right in table 6.3 and making 
a choice in each column^. 

As demonstrated in section 5.4.3 of chapter 5, when performing the Ade- 

*N.B. The Form Factor column is not applicable to the naive fits. 



Approach Form Factor Chiral Extrapolation Treatment of Lattice Lattice Spacing 

Spacing Artifacts in set from 



Adelaide 


dipole 


Cubic 


Co term has 


To 


i.e. eq.5.13 


Gaussian 


i.e. 0{Ml,s) included 


0{a) corrections 




Naive 




Quadratic 


ao term has 


a 


i.e. eq.5.14 




i.e. no 0{M^s) term 


(9(a^) corrections 





Table 6.3: The different fit types used in the global analysis. 



6.4. FITTING ANALYSIS 



125 



laide fits we have to determine a best value for tlie A parameter. We remind 
tlie reader tliat A is a lengtli scale and it is tliis parameter tliat models the 
size of the quasi-particle that we are studying. Hence it is this that controls 
the chiral physics. Again because we use a subtracted style of fit it is not 
possible to allow A to be a free parameter in the fit. We have, in our early 
works, studied fitting functions that allow A to be a free parameter in the fit 
but found that these fits were highly unstable. Instead we use the method 
outlined in section 5.4.3. This is where we manually vary the lambda pa- 
rameter then plot the x^/d.o.f for each different fitting function against A. 
Figures 6.13 and 6.14 represent this for the dipole and Gaussian form factors 
respectively. 
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Figure 6.13: A plot of x^/d.o.f against A for the dipole form factor. The 
dashed horizontal line represents increasing from its minimum value by 
unity for the tq data (i.e. it represents one standard deviation). The intercept 
of this dashed line with the curves (at A =535 and 626 MeV) is used to 
derive upper and lower bounds for the preferred A value for the dipole case. 



When the scale is set from tq we see that the x^/d.o.f as a function of 
A has a similar functional form for all fitting functions. The plots also show 
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400 500 600 700 

A [MeV] 

Figure 6.14: A plot of /d.o.f against A for the Gaussian form factor. The 
dashed horizontal line represents increasing from its minimum value by 
unity for the tq data (i.e. it represents one standard deviation). The intercept 
of this dashed line with the curves (at A =486 and 562 MeV) is used to 
derive upper and lower bounds for the preferred A value for the Gaussian 
case. 
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that the preferred value of A has a small dependency on the order of the chiral 
expansion. For the dipole case this is very small of the order of 1% and for 
the Gaussian case it is nearer 5%. Note though the change in x^/d.o.f as 
A is varied is less for the Gaussian case. We observe no dependence of the 
preferred value of A on the modelling of the lattice systematics. We also see 
that the better x^/d.o.f is given when the dipole form factor is used. As 
expected we see that the preferred value of A is dependent on the type of 
form factor used. When setting the scale from tq our preferred values of A 
are 

r \dipole ^ gQQ 

Quadratic <^ "° [MeV] 

Cubic <^ ;° . [MeV] (6.15) 

[_ A 525 . 

When the scale is set from a the plots again show a similar functional 
form for the /d.o.f for all fitting functions. We see that the preferred value 
of A again has a dependence on the order of the chiral expansion which is 
about 5%. As in the case where the scale is set from ro the preferred value of 
A has no dependence on the modelling of the lattice artifacts. When setting 
the scale from a our preferred values of A are^ 

f ^dipole _ ggg 

Quadratic <^ [MeV] 

1^ ^Gaussian 500 

f \ dipole _ t^O'i 

Cubic <^ " [MeV] (6.16) 

l_ A '^^^^'^'^ 475 

Figures 6.13 and 6.14 show that as expected the preferred value of A is 
dependent on the form factor used, but it is also dependent on the method 
used to set the scale. This discrepancy was discussed in section 5.5.1 of 
chapter 5 in terms of the mesonic data. We will address the case for the 

^N.B. We do not quote errors for the preferred values of A for the quadratic case or 
when a is used to set the scale because our final prediction for the nucleon mass will come 
from the cubic rg case. This is because the lowest x^/d.o.f is found when rg is used to set 
the scale. 
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nucleons in section 6.6. 

We see that (as with the mesonic data in chapter 5), when the scale is 
set from tq, 0{a^) errors dominate. But when the scale is set from a, 0{a) 
errors dominate. This is further evidence supporting the idea that in cases 
where the scale is set using a the dominant lattice spacing systematics will be 
0{a) and in cases where the scale is set from ro the dominant lattice spacing 
systematics will be (9(a^). As in chapter 5 we have no explanation as to why 
this should be the case. 

We use these preferred values of A to perform the 16 Adelaide fits outlined 
in table 6.3 and the 8 naive fits that are also listed in table 6.3. The results 
of these fits along with the /d.o.f for each fit are listed in tables 6.4 and 
6.5 for the cases where the scale is set from tq and a respectively. 

As with our global study of the mesonic data (sec 5.4.3) we conclude this 
section by summarising the results of the global fitting analysis (tables 6.4 
& 6.5 and figs 6.13 & 6.14). 

• Fit approach 

We see that the best /d.o.f. (indicating the best fitting procedure) is 
given by the Adelaide method which uses a dipole form factor. This has 
the best /d.o.f. in every case. This is further supporting evidence for 
the Adelaide method being a valid chiral extrapolation procedure. This 
is true for both methods of setting the scale. When a Gaussian form 
factor is employed we see that the x^/d.o.f. for the Adelaide method 
is roughly equal to that of the naive fit when the scale is set from tq. 
When the scale is set from a we see that the naive method performs 
slightly better than the Adelaide method when a Gaussian form factor 
is used. We believe that this indicates that the Gaussian form factor 
does not represent the continuum behaviour of the pions in the quasi- 
particle correctly. 

• Chiral extrapolation 

Errors in the higher order coefficients are large for the cubic fits com- 
pared to their quadratic counterparts. This said, the cubic fits always 
produce a non-zero coefficient indicating the need for a cubic term 
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Table 6.4: The results of the global fit analysis where the scale is set from tq. 
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Table 6.5: The results of the global fit analysis where the scale is set from a. 
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in m^g^. (Note that the quadratic fits were preferred in the mesonic 
case.) 

• Treatment of the lattice spacing systematics and the fit coefficients 
We note the remarkable agreement between the coefficients of the Ade- 
laide dipole and Adelaide Gaussian fits for each fitting procedure. More- 
over the coefficients of the Adelaide fits agree for the same order of the 
chiral expansion (cubic and quadratic fits). We see that the Xn coeffi- 
cients differ between the various fitting methods. We believe this is an 
indication that the lattice spacing error contained in the coefficient 
is more complicated that a simple 0{a) or 0{a'^) error. Though we 
still believe that 0{a) errors dominate when the scale is set from a and 
C(a^) errors are dominant when the scale is set from tq. 

• Setting the scale 

As seen in our study of the mesonic data (sec 5.4.3), the x^/d.o.f. is 
almost halved when the scale is set from tq compared to when the 
scale is set from a (figs 6.13 & 6.14). This, along with strong evidence 
for 0{a) lattice systematics in the a data, gives us reason to favour 
setting the scale from tq. We also note that in this study the naive 
method seems to prefer setting the scale from ro, whereas in our study 
of the mesonic data (sec 5.4.3) the naive method seemed to have no 
preference. 

Using our results from this section we select the cubic chiral extrapolation 
method with an O(a^) correction in the Oq coefficient with the scale set from 
ro to define the central values of the Adelaide and naive methods. We favour a 
dipole form factor for the Adelaide method. The spread from the other fitting 
types is used to define the error. Section 6.5 contains physical predictions 
from the nucleon mass. 

6.4.4 Analysis of the different form factors 

In this subsection we include a brief study of the two form factors that are 
employed in this chapter (dipole eq 6.4 & Gaussian eq 6.5). We do this in an 
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attempt to prove that the Adelaide method can employ different functional 
forms for the form factor u{k) and still produce similar results (within errors). 
We note though that the results from the previous subsection (sec 6.4.3) 
indicate a preference for the dipole form factor. Figure 6.15 shows a 2D plot 
of the two form factors with A set equal to the preferred values for each form 
factor (eqs 6.15 & 6.16). 
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Figure 6.15: A plot of the dipole and Gaussian form factors. Calculated at 
the preferred values of A for each form factor. 



The plot shows for the preferred values of A the functional forms of the 
dipole and Gaussian form factors are similar. It is clear though the dipole 
provides a sharper cut-off compared to the Gaussian. For a more intuitive 
view of how the form factors behave we include 3D plots for the dipole and 
Gaussian form factors (figs 6.16 & 6.17) showing how their functional form 
behaves as A and change. 

We see that while the two plots show similar behaviour for smaller values 
of A, their behaviour changes as A increases. For large values of A the dipole 
provides a far sharper cut-off compared to the Gaussian form factor. It would 



6.4. FITTING ANALYSIS 



133 




134 



CHAPTER 6. AN ANALYSIS OF THE NUCLEON MASS 




6.4. FITTING ANALYSIS 



135 



be reasonable to assume that this is why the preferred values of A are smaller 
in the Gaussian case compared to the dipole case. This provides evidence 
for the data selecting a form factor which has a functional form that most 
closely represents the continuum behaviour of the chiral physics. Since the 
data cannot alter the function which is used as the form factor it changes the 
A parameter to suit. We offer further evidence for this in the form of figures 
6.18 and 6.19. Here we see the behaviour of the self energy when the scale 
is set from tq. Figure 6.18 represents the self energy data for a dipole form 
factor with the A parameter set to our preferred value for this case which is 
594 [MeV]. Figure 6.19 represents the self energy data for a Gaussian form 
factor with the A parameter set to our preferred value for this case which is 
525 [MeV]. 
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Figure 6.18: This plot shows how the self energy behaves when a dipole form 
factor is used and the scale is set from tq. We use our preferred value of A 
in this case which is 594 [MeV]. 



Although different form factors have been used, the resulting functional 
form of the self energy appears to be very similar. Importantly the energy 
scale that the self energy covers is almost identical. We see that the data for 
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Figure 6.19: This plot shows how the self energy behaves when a Gaussian 
form factor is used and the scale is set from tq. We use our preferred value 
of A in this case which is 525 [MeV]. 
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the dipole case has a smaller spread and this may contribute to the reduction 
in the /d.o.f in our fits. 

We conclude this section by noting that the Adelaide method seems to 
prefer a value for the A parameter that causes the form factor to represent 
the chiral physics in the continuum. 

6.5 Physical predictions 

In this section we make a physical prediction for the continuum nucleon mass 
(Mat). This is done for the Adelaide method (for both types of form factor) 
and also for the naive method (eqs 6.13 & 6.14). All of our predictions in 
this section will come from our global approach (studied in section 6.4.3) 
rather than the individual approach that was employed in section 6.4.2. We 
choose to use the global method because we expect the coefficients from this 
method to be more accurate than those from the individual approach since 
these fits should be highly constrained. Here we also include a study of the 
nucleon mass prediction (Af/v) as a function of A. We obtain our continuum 
predictions by setting Mp'^ = Mp°g~'^'^^ = Mp^g^ = /i^ in equations 6.13, 
6.14 and 6.7 with /i^r being the physical pion mass which we take to be 138 
[MeV]. We also set M^^^ = M^'"'-'''' and M^^^ = M^"""'"'^ in equation 6.7. 
In doing this we see that the rj' contributions to the total self energy (eqs 

6.6 & 6.7) disappear in the continuum as required. The only remaining term 
involving and Ma is the crj^^ self energy term. We set this equal to the 
physical mass splitting of the nucleon and A which we take to be 293 [MeV] 
[33]. To calculate the self-energy terms in the continuum we use the same 
fourth order Runge-Kutta method that we employed in chapter 5. To make 
a physical prediction for each different fitting method (table 6.3) we use the 
coefficients (oq"""*, 02, (34 & ag) listed in tables 6.4 and 6.5 for the cases where 
To and a are used to set the scale respectively. We list the results for our 
physical predictions in table 6.6. We use the relevant preferred value of A for 
each Adelaide fit (taken from equations 6.15 & 6.16 for the scale set from rg 
and a respectively). 

In figures 6.20 and 6.21 we present a graphical representation of our study 
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Estimate Form 


Mn [GeV] 


Mn [GeV] 


Approach Factor 


(Scale from vq) 


(Scale from a) 


Experimental 


0.939 


0.939 



Cubic chiral extrapolation oq contains 0{a) 



Adelaide dipole 0.98411^ 0. 9501^1 

Adelaide Gaussian 0.973111 0.93811^ 
Naive - lM6t{l 0.9921^3 



Cubic chiral extrapolation oq contains (9(a^ 



11 

-12 



Adelaide dipole 0.9Q5t\l 0.934 

Adelaide Gaussian 0.9561^^ 0.923iii 
Naive - 1.023t{l 0.9741}! 



Quadratic chiral extrapolation oq contains 0{a) 

Adelaide dipole 1.0061^ 0.9741^ 

Adelaide Gaussian 0.9861| 0.9591^ 
Naive - 1.076+1 1.027+^ 



Quadratic chiral extrapolation qq contains 0{a^ 

Adelaide dipole 0.9881? 0.9581^ 

Adelaide Gaussian 0.9691? 0.9451^ 
Naive - 1.0541? 1.0081| 



Table 6.6: Estimates of Mn obtained from the global fits. Our experimental 
estimate comes from a simple average of the proton and neutron masses. 
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into the variation of M^r with the A parameter for the dipole and Gaussian 
form factors respectively. Each figure contains eight data sets corresponding 
to the four different types of Adelaide fit and the two ways to set the scale. 
For each of these plots we include an acceptable range for the A parameter 
which is represented by two vertical dashed lines. To find this range we use 
our plots of x^/d-o.f against A (figs 6.13 & 6.14). We define the range of 
acceptable values of A by increasing from its minimum by unity. This 
represents one standard deviation. For our preferred fitting method^ we 
recall from 6.15 the acceptable range for A 

jydipole ^ 594+32 
j^Gaussian ^ ^^5^11 (6.17) 



We summarise the results of this section which are outlined in table 6.6 
and figures 6.20 and 6.21 below. 

• The statistical errors in the mass estimates are typically less than 1% 
for the quadratic extrapolations and less than 2% for the cubic extrap- 
olations. 

• We see disagreement between all types of fit when different methods are 
used to set the scale. When the scale is set from tq the mass predictions 
are always higher than when the scale is set from a. We observed a 
similar effect in the case of the of the mesons (5.9), although it is less 
pronounced. 

• We also see that the mass predictions for a particular method (i.e. the 

Adelaide dipole, Adelaide Gaussian or naive method) have a variation 

in the results of between 3% and 5%, with the largest variation in 

the naive mass predictions. This disagreement suggests instability in 

the fits. We believe this is because the lattice systematics are more 

complicated than those we have uncovered. 

^We remind the reader that this is the cubic chiral extrapolation method with an 0{a^) 
correction in the ao coefficient with the scale set from ro (sec 6.4.3). 
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Figure 6.20: A plot of Mn as a function of A from the Adelaide approach 
using a dipole form factor. Recall that the best A value when the scale is set 
from ro((j) for the dipole form factor is A = 594(523) MeV. The two vertical 
dashed lines define the range of acceptable A values (535 MeV < A < 626 
MeV) obtained by increasing by unity in fig. 6.13. 
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Figure 6.21: A plot of Mn as a function of A from the Adelaide approach 
using a Gaussian form factor. Recall that the best A value when the scale 
is set from ro{a) for the Gaussian form factor is A = 525(475) MeV. The 
two vertical dashed lines define the range of acceptable A values (486 MeV 
< A < 562 MeV) obtained by increasing by unity in fig. 6.14. 
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• The Adelaide method always produces the mass prediction closest to 
the physical nucleon mass. For the cubic fits the Adelaide mass pre- 
dictions are very accurate compared to their naive counterparts, they 
are typically within two statistical errors of the experimental mass. For 
the Adelaide method to reproduce the experimental mass prediction a 
rescaling of around 3% in tq and 1% in ^/a is needed. 

• The variation of in the region of allowed values of A is very small 
for each different fit. Typically of the order of the other uncertainties. 

As with the results of the meson study (chap 5) we conclude by noting 
all of these points favour the Adelaide approach over the naive method. 
We believe that the Adelaide method should be the preferred method when 
performing chiral extrapolations and is a significant improvement over the 
naive method. To give the final value for Mj^ for both the Adelaide method 
and the naive method, we use our preferred fitting function (the cubic with 
0(0^) corrections in oq) and our preferred method for setting the scale (from 
To). For the Adelaide method we use our preferred form factor which is the 
dipole form factor. We quote an error that is based on the spread in the 
mass predictions (for the ro case only). We also (for the Adelaide method) 
include an estimate of the error associated with the A parameter which is 
taken from the vertical dashed lines in figure 6.20. 

Hence our final mass estimate for the nucleon is 

^Adelaide ^ 965(15)l«l« [MeV] (6.18) 
M^aive ^ 1023(15)lf [MeV] (6.19) 

where the first error is statistical and the second is taken from the fit pro- 
cedure. The third error in the Adelaide case is that associated with the A 
parameter. We have not considered any error that may be associated with 
the determination of tq. We see that although the Adelaide prediction has a 
slightly wider error range, at its lower limit it comes within 3 [MeV] of the 
experimental value confirming the Adelaide method to be the better chiral 
extrapolation procedure. 
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6.6 Setting the scale 

In this chapter we have studied two methods of setting the scale, from ro 
and from a. We remind the reader of the resuhs of section 5.5.1 where we 
investigated the ratio of these scales. When plotting ao-/a,.Q against we 
found that the plot had a roughly constant value of 5% above unity (fig 5.17). 
We suggested that an explanation for this could be if the ratio \^rQ = 440 
[MeV] X 0.49 [fm] is about 5% below its true value. This would explain why, 
when setting the scale from tq, the estimates of the nucleon mass (M^r) are 
always larger than those where the scale if set from a (table 6.6) for the naive 
case. For the Adelaide method the relationship between the lattice scale 
and Mjv is a highly non-linear one due to the functional form of S^ot- So we 
cannot imfere a similar relationship between mass estimates from different 
methods used to set the scale as in the naive case. We also believe that this 
non-unit ratio is the cause of the difference in preferred values of A that is 
observed when different methods are used to set the scale (figs 6.13 & 6.14). 

6.7 Conclusions 

We conclude by listing the results of our study. 

• We have shown the Adelaide method to be a valid method for chiral 
extrapolations. 

• We have applied the generalised Adelaide chiral ansatz for the nucleon 
to the "pseudo-quenched" case (i.e. when Hsca. 7^ ^^vai)- 

• We have tried to uncover unquenching effects in the data but found 
little evidence of this and have not managed to quantify them. 

• We have tried to quantify the residual 0{a) effects, but feel that we 
have not uncovered the full lattice spacing systematics. 



We have studied different fitting approaches (sees 6.4.2 & 6.4.3) and 
found that our global procedure to be the more robust method. 
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• We have demonstrated the Adelaide method can predict a preferred 
value of the A parameter with resulting errors approximately equal to 
other statistical errors in the procedure. 

• We have shown the Adelaide method can predict a preferred functional 
form. By altering the A parameter the Adelaide method causes the 
form factor to describe continuum physics as best it can. 

• We have listed 24 different predictions for the nucleon mass (sec 6.6). 
We find the Adelaide method with a dipole form factor using a cubic 
fit with C(a^) corrections in the qq coefficient to be the best fitting 
procedure. 

• We have indicated that small errors in ro and a may cause incorrect 
central values for our mass estimates (sec 6.6). 

• Finally we note that theoretically our fit procedure could be improved 
if A and the physical mass splitting between the A and nucleon were 
allowed to be free parameters in our fit procedure. This would (in 
theory) give a accurate A parameter and also allow the A mass to 
be determined. It may also be possible to use an iterative procedure 
whereby the results of a fit are used as the physical mass splitting in an 
attempt to produce a more accurate determination of the nucleon mass, 
this in turn would be used to produce a more accurate determination 
of the A mass, ad infinitum. 

As in the mesonic case we haven't modelled finite-size effects because we 
don't have enough different volumes to do this. Also although we haven't per- 
formed an infinite volume extrapolation finite volume effects are considered 
by the Adelaide method (eq 5.5) since the momentum integral is replaced by 
the appropriate kinetic sum. 



Chapter 7 
Conclusions 



This thesis describes the results of investigations into chiral extrapolation 
procedures. Our research focused on two types of extrapolation procedure, 
namely the Adelaide method and a standard polynomial extrapolation pro- 
cedure which we refer to as the naive method. These methods were used to 
produce physical mass predictions which were then compared to experimental 
results. In the mesonic case we also studied the J-parameter [21] 

To do this we simulated data which was first produced by the CP-PACS 
collaboration [22]. In the case of the Adelaide method self energy values 
corresponding to the CP-PACS data set were calculated. Our research has 
been restricted to the "pseudo-quenched" case (i.e. when Ksca 7^ ^^vai)- 

During our investigations we have 

• quantified the residual 0{a) effects in the CP-PACS data. 

• introduced a global fitting method which allows us to treat data gen- 
erated on lattices which have different lattice spacings as a single data 
set (for example sec 5.4.3). 

• demonstrated how the Adelaide method can predict a preferred value 
for the A parameter. 

• indicated that small errors in the values of tq and a might be the cause 
of the slight inaccuracy in the central values of our mass estimates. 
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• studied different form factors and shown how the Adelaide method can 
pick a preferred function. 

We now hst our final results for studies of the p-mass, the nucleon mass 
and discuss the results for the J-parameter. 

Our final result for the p-mass was obtained from the Adelaide method 
using our preferred fitting function (the quadratic fit with O(a^) corrections 
in the ao coefficient) and our preferred method for setting the scale (from 
Tq). We find 

^Adelaide ^ 779 (4) t jMeV] (7.1) 

We remind the reader that the first error is a statistical error, the second is 
associated with the choice of fitting function and the third is that which is 
related to the determination of the A parameter. The central value of our 
final estimate is just 9 [MeV] away from the experimental value (770 [MeV]) 
and is equal to the experimental value within errors. 

Our final result for the nucleon mass comes from the Adelaide method 
where we employ a dipole form factor, we use our preferred fitting func- 
tion (for the nucleon this is the cubic with O(a^) corrections in ao) and our 
preferred method for setting the scale (which is again from from tq). 

^Adelaide ^ 965(15)t«t« [MeV] (7.2) 

As with the p-mass prediction the first error is statistical, the second is 
taken from the fit procedure and the third error is that associated with the 
A parameter. The nucleon mass prediction at its lower limit it comes within 
3 [MeV] of the experimental value (939 [MeV]) 

In both cases we do not consider any error that may be associated with 
the determination of Vq. 

For the J-parameter we study J'^**^''^*^. This is because J'^^'^^rete 
easily be determined from experimental data, but J itself cannot. We remind 
the reader that table 5.9 lists values for j^^^^crete^ 'jj^ggg values of J'^*^'^'''^*^ are 
underestimates of the experimental value. This is a well known phenomena 
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and is of no surprise. 

We conclude this work by reiterating that the Adelaide method appears 
to be a valid chiral extrapolation procedure and should be favoured over 
standard polynomial fitting methods. We believe that for the foreseeable 
future the Adelaide method will prove to be a valuable tool for the extrap- 
olation of lattice data. This is because it will be many years before high 
performance computers can run lattice simulations with quark masses near 
the chiral limit. 



CHAPTER 7. CONCLUSIONS 



Bibliography 



[1] S. Sakata., 1956, Prog. Theor. Phys. 16, 686. 

[2] M. Ikeda, S. Ogawa, Y. Ohnuki., 1959, Prog. Theor. Phys. 22, 715. 

[3] M. Gell-Mann, Y. Ne'eman., 1964, The Eightfold Way. Reading: W. 
Benjamin. 

[4] K. Nishijima., 1955, Prog. Theor. Phys. 13, 285. 

[5] M. Gell-Mann., 1956, Nuovo. Cim. Supp. 4, 848. 

[6] M. Creutz., 1983, Quarks, gluons and lattices (Cambridge). 

[7] J. Smit., 2002, Introduction to Quantum Fields on a Lattice (Cam- 
bridge) . 

[8] B. Efron., 1979, Bootstrap methods: another look at the jackknife, Ann. 
Stat. 29, 1-26. 

[9] M. Gell-Mann, M. Levy., 1960, Nuovo Cim. 16, 705. 

[10] J. Gasser, M. E. Sainio, A. Svarc, 1988, Nucl. Phys. B307, 779. 

[11] S. Weinberg., 1979, Physica 96A, 327. 

[12] R.D. Young, D.B. Leinweber, A.W. Thomas., 2003, 
Prog.Part.Nucl.Phys. 50 399-417. hep-lat/0212031 

[13] B. Holstein., Chiral Perturbation Theory: A Primer, hep-ph/9510344. 

[14] B. Borasoy., 1999, Phys.Rev. D59, 054021. 

149 



150 BIBLIOGRAPHY 

[15] E. Jenkins & A. V. Manohar, Phys.Lett. B255, (1991) 558. 

[16] S. Scherer., Introduction to Chiral Perturbation Theory. 
hep-ph/0210398. 

[17] J. N. Labrenz & S. R. Sharpe., Phys.Rev. D54, (1996) 4595. 
hep-lat/9605034 

[18] C.K. Chow and S.J. Key., 1998, NucLPhys. B528, 303. 
hep-ph/9708432. 

[19] V. Bernard, N. Kaiser, Ulf-G. MeiBner., Chiral Dynamics in Nucleons 
and Nuclei, hep-ph/9501384. 

[20] Craig McNeile., Meson and Baryon Spectroscopy on a Lattice. 
hep-lat/0307027 

[21] R. Sommer., 1994, Nucl. Phys. B411, 839. hep-lat/9310022. 

[22] CP-PACS Collaboration, A. A. Khan et al, 2002, Phys.Rev. D65 
054505. hep-lat/0105015. 

[23] CP-PACS Collaboration, S. Aoki et al, 1999, Phys.Rev. D60 114508. 
hep-lat/9902018. 

[24] D.B. Leinweber et al, 2001, Phys.Rev. D64 094502. hep-lat/0104013. 

[25] C.R.Allton, V.Gimenez, L.Giusti, F.Rapuano., 1997, NucLPhys. B489 
427-452. hep-lat/9611021. 

[26] UKQCD Collaboration: C.R. Allton, et al., 2002, Phys.Rev. D65 
054502. hep-lat/0107021. 

[27] UKQCD Collaboration, P. Lacock and C. Michael., 1995, Phys. Rev. 
D52 5213. hep-lat/9506009. 

[28] R. D. Young, D. B. Leinweber, A. W. Thomas., private communication. 

[29] R. D. Young, D. B. Leinweber, A. W. Thomas, S. V. Wright., 2002, 
Phys.Rev., D66 094507. hep-lat/0205017. 



BIBLIOGRAPHY 



151 



[30] D. B. Leinweber, A. W. Thomas, K. Tsushima, S. V. Wright., 2000, 
Phys.Rev. D61 074502. hep-lat/9906027. 

[31] Malcolm Butler, Martin Savage, Roxanne Springer., 1993, Nucl.Phys. 
B399 69-88. hep-ph/9211247. 

[32] W. Armour, C. R. AUton, D. B. Leinweber, A. W. Thomas, R. D. 
Young., hep-lat/0510078 



[33] Particle Data Group., http://pdg.lbl.gov 



